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Abstract. We study the Glauber dynamics for the (2+1)d Solid-On-Sohd model above a hard 
■ wall and below a far away ceiling, on an L x L box of with zero boundary conditions, at 

D ' large inverse-temperature /3. It was shown by Bricmont, El-Mellouki and Frohlich [12] that 

pL^ , the floor constraint induces an entropic repulsion effect which lifts the surface to an average 

height H x (1//3) log L. As an essential step in understanding the effect of entropic repulsion on 
the Glauber dynamics we determine the equilibrium height H to within an additive constant: 
H = (1/4/3) log L+0{1). We then show that starting from zero initial conditions the surface rises 
' to its flnal height H through a sequence of metastable transitions between consecutive levels. 

I The time for a transition from height h = aH, a £ (0, 1), to height ft + 1 is roughly exp(cl/") for 

^ . some constant c > 0. In particular, the mixing time of the dynamics is exponentially large in L, 

' i.e. Tmix e'^^ ■ We also provide the matching upper bound Tmix ^ e'^ ^ , requiring a challenging 

^ - analysis of the statistics of height contours at low temperature and new coupling ideas and 

techniques. Finally, to emphasize the role of entropic repulsion we show that without a floor 
constraint at height zero the mixing time is no longer exponentially large in L. 



1. Introduction 



(N 
> 

. The (d-l- l)-dimensional Solid-On-Solid model is a crystal surface model whose definition goes 

\^ I back to Temperley [51] in 1952 (also known as the Onsager-Temperley sheet). Its configuration 
space on a finite box A C Z"^ with a floor (wall) at 0, a ceiling at some n"*" and zero boundary 

. conditions is the set „+ of all height functions rj on Z*^ such that A B x ^ r]x £ {0,1, . . . , n+j 

CN I whereas r]x = for all x ^ A. The probability of G ^A,n+ is given by the Gibbs distribution 

71"a(??) = ^exp ( -/3 V Ir/a; -r/yl ) , (1.1) 

^ ' where /3 > is the inverse-temperature, x ~ y denotes a nearest-neighbor bond in the lattice 

- - - Z*^ and the normalizing constant Z\ is the partition-function. 

Numerous works have studied the rich random surface phenomena, e.g. roughening, localiza- 
tion/delocalization, layering and wetting to name but a few, exhibited by the SOS model and 
some of its many variants. These include the discrete Gaussian (replacing \r]x — ^y \ by \r]x — Vyl'^ 
for the integer analogue of the Gaussian free field), restricted SOS (nearest neighbor gradients 
restricted to {0,ibl}), body centered SOS [G], etc. (for more on these flavors see e.g. [1,4,10]). 

Of special importance is SOS with d = 2, the only dimension featuring a roughening transition. 
Consider the SOS model without constraining walls (the height function r] takes values in Z). 
For d = 1, it is well known ([25,51,52]) that the SOS surface is rough (delocalized) for any /3 > 0, 
i.e., the expected height at the origin (in absolute value) diverges in the thermodynamic limit 
|A| — )• oo. However, for d ^ 3 a Peierls argument shows that the surface is rigid (localized) for 
any /3 > (see [l-]]), i.e., [ryoj is uniformly bounded in expectation. This is also the case for d = 2 
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and large enough f3 ([11,32]). That the surface is rough for d = 2 at high temperatures was 
estabhshed in seminal works of Frohlich and Spencer [29-31]. Numerical estimates for the critical 
inverse-temperature /3r where the roughening transition takes place suggest that /3r ~ 0.806. 

One of the main motivations for studying an SOS surface constrained between two wahs, both 
its statics and its dynamics, stems from its correspondence with the Ising model in the phase 
coexistence region. For concreteness, take a box of side-length L in Z,^ with minus boundary 
conditions on the bottom face and plus elsewhere. One can view the (2 -|- 1)d SOS surface 
taking values in {0, . . . , L} as the interface of the minus component incident to the bottom face, 
in which case the Hamiltonian in (1.1) agrees with that of Ising up to bubbles in the bulk. At 
low enough temperatures bubbles and interface overhangs are microscopic, thus SOS should give 
a qualitatively correct approximation of Ising (see [3,25,48]). Indeed, in line with the (2-|-1)d 
SOS picture, it is known [5] that the 3d Ising model undergoes a roughening transition at some 
/3r satisfying /3c (3) ^ ^ /3c (2) (where Pc{d) is the critical point for Ising on Z'^), yet there is 
still no rigorous proof that /3r > /3c(3) (see [1] for more details). 

When the (2 -|- 1)d SOS surface is constrained to stay above a hard wall (or floor), Bricmont, 
El-Mellouki and Frohlich [12] showed in 1986 the appearance of the entropic repulsion: for large 
enough /3, the floor pushes the SOS surface to diverge even though /3 > /Sr. More precisely, 
using Pirogov-Sinai' theory (see the review [-17]), the authors of [12] showed that the SOS surface 
on an L X L box rises, amid the penalizing zero boundary, to an average height H{L) satisfying 
(l/C/3)logL ^ H{L) ^ (C//3)logL for some absolute constant C > 0, in favor of freedom to 
create spikes downwards. 

Entropic repulsion is one of the key features of the physics of random surfaces. This phenom- 
enon has been rigorously analyzed mainly for some continuous-height variants of the SOS model 
in which the interaction potential \r]x — is replaced by a convex potential V{r]x — i]y)\ see 
e.g. [7-9, 19,55,56], see also [2] for a recent analysis of the wetting transition in the SOS model. 
As we will see below, entropic repulsion has a profound impact not only on the equilibrium shape 
of the surface but also on its time evolution under natural Markovian dynamics for the interface. 
The rigorous analysis of these dynamical effects of entropic repulsion will be the central focus 
of this work. 

The dynamics we consider is the heat bath dynamics, or Gibbs sampler, for the equilibrium 
measure tta, i.e. the discrete time Markov chain where at each step a site x G A is picked at 
random and the height r]x of the surface at x is replaced by a random variable r]'^ G {0, . . . , n"*"} 
distributed according to the conditional probability T^hi'lily^ U This defines a Markov 

chain with state space „+, reversible with respect to tta, commonly referred to as the Glauber 
dynamics. As explained below, our results apply equally well to other standard choices of 
reversible Markov chains, such as e.g. the Metropolis chain where only moves of the type r/^ = 
% ± 1 are allowed. 

The mixing time Tmk is defined as the number of steps needed to reach approximate station- 
arity with respect to total variation distance, see Section 2 for definitions. 

The main result of this paper is that the mixing of Glauber dynamics for the (2-|-1)d SOS is 
exponentially slow, due to the nature of the entropic repulsion effect. 

Theorem 1. For any sufficiently large inverse-temperature f3 there is some c{/3) > such that 
the following holds for all L G N. The mixing time T^ix of the Glauber dynamics of the (2 -|- 1)d 
SOS model on A = {I, . . . ,L}^ with zero boundary conditions, floor at zero and ceiling at 
with log L ^ ^ L satisfies 



(1.2) 
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Figure 1. Illustration of the series of metastable states in the surface evolution. 
The dynamics waits time e^^^pi'^l^f^) until the formation of a macroscopic droplet 
(marked in red) which eventually raises the average height from h — 1 to h. 

The exponentially large mixing time in (1.2) is in striking contrast with the rapid mixing 
displayed by Glauber dynamics of the (1 + 1)D SOS model [15,44]. When d = 1 it is known 
that the main driving effect is a mean-curvature motion which induces a diffusive relaxation to 
equilibrium, with Tmk of order up to poly(logL) corrections. As we will see, in (2 + 1)d 
instead the main mechanism behind equilibration is a series of metastable transitions through 
an increasing series of effective energy barriers caused by the entropic repulsion. This is also in 
contrast with the behavior of related interface models with continuous heights as e.g. in [20,24]. 

1.1. Metastability and entropic repulsion. Consider the evolution of an initially flat surface 
at height zero. We shall give a rough description of how it rises to the final height H{L) through 
a series of metastable states indexed by /i ^ 0. Roughly speaking the surface in state with 
label h is approximately flat at height h with rare up or downward spikes. Of course downward 
spikes cannot be longer than h because of the hard wall. If /i < H[L) then the surface has an 
advantage to rise to the next level h + 1. This is due to the gain in entropy, measured by the 
possibility of having downward spikes of length /i + 1, beating the energy loss from the zero 
boundary conditions. 

The mechanism for jumping to the next level should then be very similar to that occurring 
in the 2d Ising model at low temperature with a small external field opposite to the boundary 
conditions (see [5-3,54]). Specifically, via a large deviation the surface at height h creates a large 
enough droplet of sites at height /i + 1 which afterwards expands to cover most of the available 
area. The energy /entropy balance of any such droplet is roug hlyi of order /3|7| - e-^^('^+iU(7) 
where I7I and ^(7) are the boundary length and area respectively and the effective field e"^'^^''"''^-* 
represents the probability of a 1 x 1 x (h + 1) isolated downward spike. Simple considerations 
suggest then that the critical length of a droplet should be proportional to e^^^^'^'^^ Finally the 
well established metastability theory for the 2d Ising model indicates that the activation time 
for such a critical droplet should be exponential in the critical length^ (i.e. a double exponential 
in h) as seen in Figure 1. 

-'^Here we are neglecting finer results taking into account the surface tension and the associated Wulff theory; 
the basic conclusions of this reasoning are nevertheless still valid. 

^At the early stages of the process when h is quite small the activation time has important corrections to this 
guess due to the many locations in the L x L box where the droplet can appear. However, as soon as h becomes 
of order log log L these entropic corrections become negligible. 
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Figure 2. Glauber dynamics for SOS on a 64 x 64 square lattice at /3 = 0.9 from 
an initial state r] = 10. Surface gradually falls towards level H = 1. Snapshots 
at t = 10 (top left), t = 100, t = 1000, t = 10000 (bottom right) in cont. time. 



Of course, in order to establish, even partially, the above picture and to prove the asymptotic 
of log(TMix) as per (1.2) it is imperative to estimate the final equilibrium height of the surface 
H{L) to within an additive 0{\). In Section 3 (Theorem 3.1), we improve the estimates of [12] 
to show that in fact the typical height of the surface at equilibrium is H{L) + 0(1), where 

H{L)=[^\ogL\. (1.3) 

The aforementioned picture of the evolution of the SOS surface through a series of metastable 
states is quantified by the following result. 

Theorem 2. For any sufficiently large inverse-temperature f3 there is some c{f3) > such that 
the following holds. Let {r]{t))t ^ o be the Glauber dynamics for the SOS model on A = {1, ... , L}^ 
with zero boundary conditions, floor at zero and ceiling at with logL ^ n"*" ^ L, started from 
the all-zero initial state. Fix a € (0, 1) and let Ta = min{t : r]{t) € Vta} where 

^a = {v^ ^A,n+ ■■ #{x : Vx ^ aH{L)} > ^|A|} . (1.4) 

Then lim/,^oo 7rA(^a) = 1 and yet 

limPfe^^"^T, ^e^V-)^''^!. (1.5) 



In fact, we prove this with the constant ^ in (1.4) replaced by 1— e(/3) where lim^^oo ^{P) = 0. 



_9_ 

10 

Moreover, the statement of the above theorem remains valid when a = a{L) — )• 1 as long as the 
target level h = aH{L) satisfies h ^ H{L) — c for some sufficiently large c(/3) > 0. 

Remark. A natural conjecture in light of Theorem 2 is that there exists a constant A such that 
the distribution of Ta x e~^^'^ converges as L —t- oo to an exponential random variable. 

We wish to emphasize that, as will emerge from the proof, the exponential slowdown of 
equilibration is a coupled eff'ect of entropic repulsion and of the rigidity of the interface. In 
particular, the following rough upper bound shows that the situation is very much different 
when the fioor constraint is absent (yet the ceiling constraint remains unchanged): 

Theorem 3. Consider the (2+l)'D SOS setting as in Theorem 1 with the exception that the 
surface heights belong to the symmetric interval [— n~^,n^]. Then Tmk ^ exp(o(L)). 
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Specifically, our proof gives the estimate T^ix ^ exp(L2+°(i)). No effort was made to improve 
the exponent ^ as we would expect the true mixing behavior to be polynomial in L. We further 
expect that in the presence of a floor yet for /3 < /3r the mixing time will have a different scaling 
with the side- length L. 

It is useful to compare our results with those of [17], where the Glauber dynamics for the 
(2 + 1)d SOS above a hard wall, at low temperature and in the presence of a weak attracting 
(towards the wall) external field was analyzed in details. There it was proved that certain 
critical values of the external field induce exponentially slow mixing while for all other values 
the dynamics is rapidly mixing. Although the slow mixing proved in [17] is similar to the one 
appearing in (1.2), the physical phenomenon behind it is very different. When an external field 
is present, a critical value of it results in two possible and roughly equally likely heights for the 
surface. In this case slow mixing arises because of the presence of a typical bottleneck in the 
phase space related to the bi-modal structure of the equilibrium distribution. In the setting of 
Theorems 1 and 2 instead, there is in general no bi-modal structure of the Gibbs measure and 
the slow mixing takes place because of a multi-valley structure of the effective energy landscape 
induced by the entropic repulsion which produces a whole family of bottlenecks. 

1.2. Methods. We turn to a description of the main techniques involved in the proof of the 
main theorems. Our results can be naturally divided into three families: equilibrium estimates, 
lower bounds on equilibration times, and upper bounds on equilibration times. 

Equilibrium Estimates. Our proof begins by deriving estimates for the equilibrium distribution 
which are crucial to the understanding of the dynamics (as discussed in Section 1.1) and of 
independent interest. Over most of the surface, the height is concentrated around H{L) as 
defined in (1.3) with typical fluctuations of constant size. Achieving estimates with a precision 
level of an additive 0(1) turns out to be essential for establishing the order of the mixing time 
exponent: indeed, analogous estimates up to some additive g{L) tending to oo with L would set 
off this exponent by a factor of e^^^\ 

The main techniques deployed for this part are a range of Peierls-type estimates for what we 
refer to as /i-contours, defined as closed dual circuits with values at least h on the sites along 
their interior boundary and at most h — 1 along their exterior boundary. In the simpler setting 
of no fioor or ceiling (i.e., the sites are free to take all values in Z as their heights), the 
which decreases all sites inside an /i-contour 7 by 1 is bijective and increases the Hamiltonian 
by I7I, the length of the contour. Hence, the probability of a given /i-contour in this setting is 
bounded by exp(— /3|7|). Iterating estimates of this form allows us to bound the deviations of 
the sites with the correct asymptotic in the setup of having no walls. 

The presence of a fioor renders this basic Peierls argument invalid since the map may 
leave sites in the interior with negative values. Rather than a technicality, this in fact lies at 
the heart of the entropic repulsion effect. We resort to estimating the probability that a given 
/i-contour has a strictly positive interior, a quantity directly involving its area. By analyzing 
an isoperimetric tradeoff between the contour's area and perimeter we show that large contours 
above height H(L) are unlikely, which in turn implies 0(1) typical fiuctuations above this level. 
For a lower bound on the typical height of the surface we show that if too many sites are below 
H{L) — k then the loss in energy due to raising the entire surface by 1 is more than compensated 
by the increased entropy from the freedom to create downward spikes reaching 0. Put together, 
these estimates guarantee that the height of most sites is within a constant of H{L). 

Equilibration times: lower bounds. Fix h = aH{L) — 1 with a G (0, 1) and consider the restricted 
ensemble obtained by conditioning the equilibrium measure on the event A that all /i-contours 
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have area smaller than SL?"", for some small (5 > 0. Our equilibrium estimates imply that in this 
restricted ensemble 

(i) each /i-contour is actually very small (e.g. with area less than log(L)^), with very high 
probability; 

(ii) the probability of the boundary of A is 0(exp(— cL"^)); 

(iii) the probability of having a large density of heights at least h + 1 = aH{L) is 0(exp(— cL")). 

In some sense (i) , (ii) and (iii) above establish a bottleneck the Markov chain must pass through 
and thus provide the sought lower bound of exp(c(/3)L'^) on the typical value of the hitting 
time Ta in Theorem 2 when the initial state is the all zero configuration. In fact, the initially 
flat configuration can be replaced by monotonicity by the restricted ensemble described above. 
Then, in order for Ta to be smaller than T, either the dynamics has gone through the boundary 
of A before T or the event described in (iii) occurred without leaving A. Either way an event 
with 0(exp(— cL")) probability occurred and the minimal time to see it must be proportional 
to the inverse of its probability. 

Equilibration times: upper hounds. By the monotonicity of the system it is enough to consider 
the chain starting from the maximum and minimum configurations. The natural approach is 
to apply the well-known canonical paths method (see [21,22,35,50] for various flavors of the 
method). As the cut- width of the cube is L^, the most naive application of this approach would 
give a bound of exp(0(L^)). A better bound can be shown by considering the problem with 
maximum height n"*" = log L. In this case the cut-width is of order L log L yielding a mixing time 
upper bound of exp(0(L log L)). Since the height fluctuations are logarithmic, we can iterate 
this analysis using monotonicity and censoring to get a bound of exp(0(L log L)) for the original 
model with = L, vs. our lower bound of exp(cL). However, removing the logL factor that 
separates these exponents entails a signiflcant amount of extra work. 

The basic structure of the proof is to flrst establish a burn-in phase where we show that, 
starting from the maximal and minimal configurations, the process reaches a "good" set featuring 
small deviations from the equilibrium level H(L). From there we establish a modified canonical 
paths estimate (Theorem 2.4), showing that it is enough to establish a reasonable probability of 
hitting the good set from any starting location together with a good canonical paths estimate 
restricted to this set. This new tool, which we believe is of interest on its own right, is described 
in detail in Section 2.3 and proved in a general context in Section 5. 

Showing that the surface falls down from the ceiling (the maximum height) to H{L), as 
depicted in Figure 2, ought to have been the easier part of the burn-in argument since high 
above the floor there is no entropic repulsion effect. Unfortunately a number of major technical 
challenges must be overcome. 

First, the effect of the entropic repulsion is still apparent for the estimates we require when 
the surface is fairly close to H{L). To overcome this we add a small external field to the model, 
thereby modifying the mixing time by a factor of at most exp(0(L)) (which is large but still of 
the same order as our designated upper bound) and tilting the measure to remove these entropic 
repulsion effects. Second, while our main equilibrium estimates were proved using Peierls-type 
estimates, for the burn-in we require some of the cluster expansion machinery of [23] which we 
extend to the SOS framework. This involves a number of challenges including showing that 
the contours we consider do not interact significantly with the boundary conditions, a highly 
non-trivial fact. Implementing this scheme is the biggest challenge of the paper and we provide 
extensive notes for the reader in these sections to explain the rather technical proofs. 

Finally, the fact that the surface rises from the fioor (the all zero initial condition) to the 
vicinity of the equilibrium height H{L) in time exp(0(L)) is proved via an unusual inductive 
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scheme. Unlike other multi-scale inductive schemes, somewhat surprisingly the one used here 
does not incur any penalizing factor on the upper bound. We first prove weaker bounds on the 
mixing time and use these estimates to show that a smaller box of side- length L / log L mixes by 
time exp(0(L)). By monotonicity, we can use this to bound the distance from the equilibrium 
height of the surface in the original box by H(L) — H{L/\ogL). By using this height estimate 
along with our canonical paths result we get improved bounds on the mixing time. This in turn 
allows us to take larger sub-boxes and iteratively achieve better and better estimates on the 
distance to H[L). After sufficiently many iterations we show that the surface reaches height 
H{L) — 0(1) in time exp(0(L)) and thereafter the canonical paths estimate completes the proof. 

1.3. Related open problems. 

Tilted walls. An interesting and to our knowledge widely open problem concerns the SOS model 
with a non-horizontal hard wall, i.e. when the constraint rjx ^ is replaced by r]x ^ 4>x, where 
0" denotes the discrete approximation of the plane orthogonal to the unit vector n, and n is 
assumed to have all components different from zero. The equilibrium fluctuations for (3 = +oo 
can be analyzed via their representation through dimer coverings [3(i] and the variance of the 
surface height in the middle of the box can be shown to be 0(log L); see [14, Section 5] for a proof. 
Moreover, at /3 = +oo, as far as the dynamics is concerned, it has been proved [15] that the 
mixing time is of order up to polylog(-L) corrections and that the relaxation process is driven 
by mean curvature motion. The case (5 < +oo, however, remains open both for equilibrium 
fluctuations and for mixing time bounds. 

Mixing time for Ising model. In view of the natural connection with the Ising model, the study 
of Glauber dynamics for the SOS can also shed some light on a, still open, central problem 
in the theory of stochastic Ising models: its mixing time under an all-plus boundary in the 
phase coexistence region. The long-standing conjecture is that the mixing time of Glauber 
dynamics for the Ising model on a box of side-length L with all-plus boundary should be at 
most polynomial in L at any temperature. More precisely, the convergence to equilibrium should 
be driven by a mean-curvature motion of the interface of the minus droplet in accordance with 
Lifshitz's law [JO]. For instance, the mixing time of Glauber dynamics for Ising on an L x L 
square lattice is conjectured [26] to be of order in continuous time. This was confirmed at 
zero temperature [18,27,38] and near-zero temperatures [14], yet the best known upper bound 
for finite f3 > /3c remains quite far, a quasi-polynomial bound of L'^(^°^^^ due to [41]. The 
understanding of 3d Ising is far more limited: while at zero temperature bounds of L'^~^°^^^ were 
recently proven in [14], no sub-exponential mixing bounds are known at any finite /3 > /3c. 

2. Definitions and tools 

2.1. Glauber dynamics for Solid-On-Solid. Let U and □ denote the minimal and maximal 
configurations in Oy^n+j i-^- = and fix = n'^ for every x G A. Given a finite connected 
subset A C Z^, let (9A denote its external boundary, i.e. the set of sites in A"^ which are at 
distance 1 from A. To extend the SOS definition to arbitrary boundary conditions (b.c.) given 
by ^ : ^ Z, define the SOS Hamiltonian with b.c. ^ to be 

^i(^)-^ E \'nx-Vy\+ Yl \^--^y\- (2-1) 

x,y&A xeA, yedA 

\x-y\=l \x-y\=l 
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Given /3 > and n"*", the Gibbs measure vrf on ^\n+ with b.c. ^ is defined as 



4(^) = 4exp[-/3?^i(r?)j. (2.2) 

Notation 2.1. In the sequel when the b.c. = n £ Z we will use the abbreviated form ir^. We 
will occasionally drop the subscript A and superscript ^ from the notation of vr^ when there is 
no risk of confusion. Moreover we will need to address the following variants of ir^: 

(i) the measure t:^ of SOS without walls (no floor and no ceiling) and with b.c. at height n; 

(ii) the measure 11^ corresponding to vr^ with n"*" = +oo (no ceiling); 

(Hi) starting from Section 6 the measures tt^-^ (and its analog II^''^ with no ceiling) corresponding 
to the SOS Hamiltonian with an additional external field of the form 'Yliy^K fiVy) with 
\f\oo = 0{e~'^^) for some fixed constant c (see e.g. (Q.l)). 

The dynamics under consideration is a discrete-time Markov chain (f?(t))t=o,i,...; defined as 
follows. To construct r]{t + 1) given rj{t), 

• pick a site x £ A uniformly at random; 

• sample a new value for r]x{t + 1) from the equilibrium measure vr^ conditioned on the 
current heights at the neighboring sites, that is, r/(t+ 1) ~ 7r^(r/ G • | rjy = riy{t) Vy / x). 

The law of the process with initial condition Q is denoted by P^, the configuration at time t is 
vi^{t) and its law is When there is no need to emphasize the initial condition we simply write 
r/(i) for the configuration at time t. It is well-known that this Markov chain is reversible w.r.t. 
the invariant measure vr^. 

The mixing time T^lx is defined to be the time the process takes to converge to equilibrium 
in total variation distance, that is, 

Tmix = inf i t > : max - tt^W > (2.3) 

where ||// — denotes the total variation distance between two measures ^^v. It is well-known 
(e.g., [39, Section 4.5]) that the total variation distance from equilibrium decays exponentially 
with rate Tmix, namely 

max W^xl - TT^II ^ e-LV^MixJ _ (^2.4) 

The relaxation time Trel is the inverse of the spectral gap of the transition kernel of the chain. 
The spectral gap, denoted by gap, has the following variational characterization: 

where P is the transition kernel of the chain, / is the identity matrix and the infimum is over 
all non-constant functions /. The following standard inequality (see, e.g., [39, Section 12.2] 
and [49]) relates the mixing time and the relaxation time: 

T'rel - 1 =^ Tmix Trel log(2e/7rmin) (2.6) 

with TTmin := min^gQ^ ^+ T^^^ivi)- By definition, in the SOS model 1^2^^^+ 1 = ('^^ + 1)'^' 
'''"min ^ exp(— 4/3|A|?i+)/|r2^^^+|, thus for large enough n+ 

Trel - 1 5^ Tmix ^ 5/3|A|n+ Tree • (2.7) 
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From now on we refer to the Markov chain defined above as the Glauber dynamics. One can 
use standard comparison estimates to obtain equivalent versions of our main results for other 
standard choices of Markov chains that are reversible w.r.t. the SOS Gibbs measures, such as e.g. 
the Metropolis chain with ±1 updates. Indeed, since the heights are confined within an interval 
of size 0{L) it is not hard to see that the the ratio between the different mixing times is at most 
polynomial in L. We refer to e.g. [14, Section 6] for a detailed argument in this direction. 

2.2. Monotonicity. Our dynamics is monotone (or attractive) in the following sense. One 
equips the configuration space with the natural partial order such that a ^ rj li ax ^ r]^ for 
every x £ A. It is possible to couple on the same probability space the evolutions corresponding 
to every possible initial condition ( and boundary condition ^ in such a way that if ^ ^ ^' and 
C ^ C' then ri<{t,C) ^ V^\t,C') for every t. Here, we indicated explicitly the dependence on the 
boundary conditions but we will not do so in the following. The law of the global monotone 
coupling is denoted P. 

A first consequence of monotonicity is that the FKG inequalities [2S] hold: if / and g are 
two increasing (w.r.t. the above partial ordering) functions, then 'iT^{fg) ^ '^\{f)'^\{g) and the 
same holds for the measure vr^ without the floor/ceiling. 

Monotonicity also implies the following standard fact (cf. e.g. the proof of [ , Eq. (2.10)]): 
for every initial condition rj and boundary condition ^, 

ll/xl' - 411 ^ 2n+|A| maxdl^^ - ||/.^ - vri||). (2.8) 

Another consequence of monotonicity is the so-called Peres- Winkler censoring inequality. Take 
integers = to < ti < ■ ■ ■ < tf^ = T , a sequence of Vi C A and ^ Oj ^ 6/ ^ n+, i ^ A:. Consider 
the following modified dynamics (??(t))o ^ t ^ T- To construct f]{t + 1) given ri{t), 

• pick a site x £ A uniformly at random; 

• at time t with < t ^ ti do as follows: 

— if X ^ Vi or ii X £ Vi and fjx{t) ^ {oj, . . . ,bi} then do nothing; 

— if X £ Vi and ^ flxit) ^ bi then replace its value with a new value f]x{t + 1) in 
{oj, . . . , 6j} with probability proportional to the stationary measure conditioned on 
the value of the neighboring columns, 

fi{t + 1) ~ 7r^(r/ G • I G {a^, . . . , fej} , % = %(i) Vy / x) . 

Call flf the law at time t when the initial distribution is v. The following then holds: 

Theorem 2.2 (special case of [4G, Theorem 1.1]). If the initial distribution v is such that 
^{v)/'^a{v) o.'n increasing (resp. decreasing) function, then jl^ (rj) / 7r^^{r]) is also increasing 
(resp. decreasing) for t and ^ jl^ (resp. jl^ ^ fi^). In addition, 

M-4\\ ^ (2-9) 

2.3. An improved path argument. Geometric techniques can prove very effective in getting 
upper bounds on the relaxation time and therefore on the mixing time of a Markov chain [21, 
22,35,50] (see also [39, Section 13.5]). Let us recall the basic principle. 

Let (^(i))f=o,i,... be a discrete-time reversible Markov chain on a finite state space il, with 
invariant measure vr. For a,b £ such that the one-step transition probability p{a, b) from a 
to b is non-zero, set Q{a,b) = ■K{a)p{a,b) = Q{b,a). For each couple (c, d) E 0^ fix a path 
7(c, (i) = (xi, . . . ,Xn) in 17 with xi = c,Xn = d and p(xj,Xj_|_i) / and let [7(0, (i)| := n. Then 
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the relaxation time of the Markov chain is bounded as 

T'rel ^ , . , n/ h\ l7(??,??')k(??)vr(r/'). (2.10) 

{a,by.Q{a,b)j^O Q{a,b) ^ 

(a,fe)67(»y,r)') 

Here, (a, 5) G 7(??>V) means that if 7(7/, r/') = (xi, . . . , Xn) then there exists i such that a = 
Xi,b = Xi+i. The proof is simply an application of the Cauchy-Schwarz inequality; see e.g. [49]. 
An application of this principle gives the following: 

Proposition 2.3. For the SOS dynamics in the A = {1, . . . , L} x {1, . . . , m},m ^ L, with floor 
at height zero, ceiling at and b.c. ^, one has for some c = c(/3) 

Trel ^ cL^m^n+ exp(7/3mn+) (2.11) 

and, thanks to (2.7), T^ix = exp(0(/3Ln"'")) if L = m. 

That (2.11) easily follows from (2.10) was observed in [13] in the case of the Glauber dynamics 
of the Ising model (in this case one refers to the paths 7(r/,7/') "canonical paths"). For SOS 
the proof is very similar and is given for completeness in Section 5.1. 

However this upper bound is too rough for our purposes since we have ^ log L while we 
wish to get a mixing time upper bound which is exponential in L. Therefore a significant part 
of the present work is devoted to getting rid of the non-physical factor n"*" in the argument of 
the exponential in the r.h.s. of (2.11). Although this task may appear to be mainly of technical 
nature it actually requires a much deeper understanding of the actual behavior of the dynamics 
compared to that provided by canonical paths, and the support of new ideas. 

One of the key ingredients we use is the following improved version of (2.10), which we believe 
can be interesting in a more general context: 



Theorem 2.4. Let G C and assume that, for some T > and for every initial condition x, 
¥^{X{T) € G) ^ a with ¥^ denoting the law of the chain starting at x. Assume further that for 
every r],r]' in G there exists a path 7(77,77') as above which stays in G and let 

W{G):= max — i— V |7(r?, r?')k(??)vr(r?'). (2.12) 



a,beG Q{a,b) 

Q(a,fe)^0 ^V6G: 

(a,b)&'y{r],rj') 



Then, 



1 6 / r2 W(G)\ 
gap-i ^ - + (2.13) 

a \Pmm a ) 

with pmin := min{p((T, o"') > : o", cr' G f]} 

This is clearly an improvement provided that a is bounded away from zero, that W{G^ <C 
M^(Q) and that T is not too large (in simple words, we need that with non-zero probability the 
chain enters "quickly" the good set G where canonical paths work well). 

In our SOS application, roughly speaking, we will choose G to be the set of configurations such 
that |Ai|~^ X^xeA Vl^ ~ upper bounded by a constant. We will see that, irrespective of 

the starting configuration, at time T = exp(0(L)) the dynamics is in G with probability at least 
i. On the other hand, a minor modification of Proposition 2.3 will give W{G) = exp(0(/3L)). 
Then, Theorem 2.4 allows us to improve the mixing time upper bound to Tmk = exp(0(/3L)). 
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3. Equilibrium results 

Theorem 3.1. Let A C be a box of side-length L and let /5 ^ 1. Set H = log LJ . There 
exist some absolute constants C,K > (with K integer) such that for any integer k ^ K, 

ttI (^#{v ■.i]^i^H-k}> 6-2/5^X2) ^ exp (-c^'^l) , (3.1) 
^0 (^#{v ■.i]^^H + k}> e-'^^^L^^ ^ exp [-Ce-'^^^L{l A e-^^'^Llog-^ L)) . (3.2) 

(Notice that the bound on downward fluctuations improves with the size of the deviation 
whereas the bound on upward fluctuations deteriorates with the distance.) 

Recah that vr^ has a floor at and a ceihng at height logL ^ n'^ ^ L (together with zero 
boundary conditions). It wih be convenient throughout this section to work in the setting of a 
floor at but no ceihng, where the corresponding measure 11^ is asymptotically equal to vr^: 

Lemma 3.2. There is an absolute constant c > such that for any (3^1 and any subset of 
configurations A C {0, . . . , n'^}^ , 

Ill{A) ^ ttI{A) ^ (l + cL^e-^/^"^) n^(^) . 

The above lemma, which will be proved further on in this section, entitles us to derive results 
on vr^ from 11^ at an asymptotically negligible cost. 

The following notion of a contour and that of an h-contour, a level line at height h, play a 
crucial role in our proofs. 

Definition 3.3. We let 1?* be the dual lattice of 1? and we call a bond any segment joining 
two neighboring sites in 7?'* . Two sites x,y in 1? are said to be separated by a bond e if their 
distance (in M?) from e is ^. A pair of orthogonal bonds which meet in a site x* G Z^* is said to 
be a linked pair of bonds if both bonds are on the same side of the forty-five degrees line across 
X* . A geometric contour (for short a contour in the sequel) is a sequence eo,...,en of bonds 
such that: 

(1) Ci Cj for i 7^ j, except for i = and j = n where cq = Cn- 

(2) for every i, Ci and ej+i have a common vertex in 1?* 

(3) i/ ej, ej+i, ej, ej+i intersect at some x* S , then ej,ei+i and ej,ej+i are linked pairs 
of bonds. 

We denote the length of a contour 7 &y I7I, its interior (the sites in 1? it surrounds) by Ay and 
its interior area (the number of such sites) by |A-y|. Moreover we let be the set of sites in 
T? such that either their distance (in M."^) from j is ^, or their distance from the set of vertices 
in 1?* where two non-linked bonds of j meet equals l/\/2- Finally we let Aj^ = n and 
A-=A,\A+. 

Definition 3.4. Given a contour 7 we say that 7 is an /i-contour for the configuration rj if 

We will say that ^ is a contour for the configuration r] if there exists h such that ^ is a h-contour 
for 1]. Finally 'T^-y^h will denote the event that 7 is an h-contour. 

To illustrate the above deflnitions with a simple example, consider the elementary contour 
given by the square of side 1 surrounding a site x S Z^. In this case, 7 is an /i-contour iS rj^ ^ h 
and r]y ^ h — 1 iov all y £ {x ± ei, x ± e2, X + ei + e2, X — ei — 62}. In general, A+ (resp. A;^) is 
the set of X G A^ (resp. x £ A^) either at distance 1 from A^^ (resp. A^) or at distance y/2 from 
a vertex y £ A^ (resp. y € A^) in the south-west or north-east direction. 
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Remark 3.5. As the reader may have noticed the definition of an h-contour is asymmetric in 
the sense that we require the minimal height of the surface at the inner boundary of j, , to 
he larger than the maximum height at the external boundary. In a sense this definition covers 
upward fluctuations of the surface. Of course one could provide the reverse definition covering 
downward fluctuations. In the sequel the latter is not really needed thanks to monotonicity and 
symmetry arguments. We also observe that, contrary to what happens in e.g. Ising models, a 
geometric contour 7 could be at the same time a h-contour and a h' -contour with h ^ h' . More 
generally two geometric contours 7, 7' could be contours for the same surface with different height 
parameters even z/7 n 7' 7^ (but one of them must be contained in the other). 

The following estimates play a key role in the proof of Theorem 3.1. 

Proposition 3.6. There exists an absolute constant Cq > such that for all P ^ 1 and h ^ 1, 

ttI {%^h) ^ exp (-/3|7| + ColA^le-^'^'*) . (3.3) 
Moreover, for any family of h-contours {(7s, /is)}s6,y such that for all i^ 1 

U ^7. ^ U 

hs=i+l hs=i 

and A^^ n A^^, = when hg = h^i , s ^ s' ,we have 

<( n '^^sAs) ^ exp ( ^ (-/3|7s| + ColA.Je-^'^^^) ) . (3.4) 

As a step towards the proof of the above proposition, we consider the setting of no floor and 
no ceiling, where the picture is simpler as there is no entropic repulsion. 

Lemma 3.7. For any h-contour 7 in any domain A with any boundary condition ^ we have 

T:i{%,h) ^ exp(-/3|7|). 
Moreover, if h' < h and 7, 7' are contours with C Ay then 

4 I %',h') ^ exp(-/3|7|) . (3.5) 

Proof. Define the map T = : ^ by 

N ivv-l ve A^, 
{Tv)v = < ^, (3.6) 

I r]y otherwise. 

If r] has an /i-contour at 7 then the difference along every edge in crossing 7 decreases by 1 

Since T is a bijection it follows that 

Equation (3.5) follows from the same argument by noting that if 77 G ^-y^h H "^y^/i' then T^r/ 
remains in ^-y'^h'- This completes the proof. □ 

Remark 3.8. In the context of considering the interior of an h-contour 7 for possibly nested 
contours (such as the ones featured in Eq. (3.5)j, a useful observation is that 
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for any boundary condition ^ that is at most h all along A~ . This follows from the fact that 
conditioning on any fixed ^ h would contribute an equal pre-f actor to all configurations thanks 
to having r/f^+ ^ h, and as this includes all ^'s with r]\ ^ h — 1 this further includes '^■y^h- 

Moreover, the same holds when conditioning on ^-y^h ^ E (instead of just ^-^^h) for an arbitrary 
event E which is only a function of the configuration on (A^)"^. 

(Note that the above remark similarly applies to 11 and tt by the same argument.) 
A Peierls-type argument will transform the above lemma into the following bound on upward 
(downward) fluctuations in the no floor, no ceiling setting. 

Proposition 3.9. There exists an absolute constant c > such that for any /3 ^ 1, domain A, 
site v G A and height h ^ 0, 



2 



Proof. Define the map 5 : Z"^ — > by {Srj)^ = rju for u ^ v and 

{Sri)y 



\r]^ + h ^ , 
\^r]y - h r]y <0. 

Observe that KS't])^, | ^ h and that since S changes the Hamiltonian by at most 4/i, 

rcliSr^) -^e-'^Hlir^) . 
Moreover, as S is injective, summing over r/ we have that 

4 (1^.1 ^ /^) = E -A (^^) ^ E -A (^) = 

Since by symmetry {rj^ ^ /i) = vr^ (ry„ ^ — h) the lower bound follows. 
To get the upper bound define a set of nested contours surrounding v as 

A{h, -y) = { (71, . . . , 7/i) : D G A^^ and A-y.^-^ C A^. for ah 1 ^ i ^ /i - l} 

and observe that, if rj is such that ry„ ^ h, then necessarily there exists (71, . . . ,7h) £ A{h,v) 

such that r/Gfli^i^ h'^Tni- 

Applying Lemma 3.7 iteratively (while bearing Remark 3.8 in mind) we now obtain that for 

every (71, . . . ,7h) S A{h,v), 

-i(ni,.,.^...) ^e-^Et.W. (3.7) 

Simple counting gives that the number of contours of length n starting from a vertex is at most 
Rn, the number of self avoiding walks of length n. If such a path surrounds v then it must cross 
the horizontal line containing v to its right within distance n so the number of 7 with I7I = n 
and G A^ is at most (with room to spare). Hence 

00 

7 : iigA^ n=8 
|A^|>2 

1/n 

which is uniformly bounded in /3 for any /3 ^ 1 since the connective constant ^2 = hnin-s^oo Rn 
is known to satisfy ^2 < 2.68 < e. Hence, for some large enough M, independent of /3, 

^ e-P\^\^e-^P. (3.8) 

7 : 1'SA^ 
|A^|>M 
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Now define a collection of nested contours of area at least 2 and at most M as 

AM{h,v) = {(71,. . . ,7h) G A{h,v) : 2 ^ |A^.| ^ M for all 1 i /i} . 
We note that 

/ /i + M — 2\ „,mA/-1/i , 7i,r\A/-2 



\AMih,v)\^\AMil,v)r~'(^ j<\AMil,vT-'ih + Mr'-' (3.9) 

since, examining the way |A^. | decreases, there are at most M — 2 transitions of |A^. | < |A^._-^| 
and in each case the number of possible ji is at most |^m(1) ^^)| with much room to spare. 

For any (71, . . . ,7/^) S A{h, v) we can find ^ k ^ I ^ h such that |A^. | > M for 1 ^ i ^ k, 
that (7fc+i, ■■■,70 £ AAiil — k,v) and |A^J = 1 iov I < i ^ h. Then 

(■yi,...,"/h)£Aih,v} Os:k<:i<:h (71,. ..,7^) [■yk+i,-;ji)eAMil-k,v) 



^ A: ^ / h 

^ e-^l"' Y |^A./(l,t;)|*'-^(/ -k + Mf'-^e-^P' 



where the first equality holds since |7i| = 4 when |A^. | = 1, the inequality in the second line 
is by equation (3.8) and the fact that every contour with |A^| ^ 2 has I7I ^ 6, and where the 
transition in the third line is by (3.9). Combining with equation (3.7) completes the proof. □ 

Proposition 3.9 allows us to readily infer Lemma 3.2: 

Proof of Lemma 3.2. One has 

vri(A) = n^(^)^ 

where denotes the partition function corresponding to H^. The fact that n^(^) ^ 7r^(^) 
follows immediately from ^ Z^. To show that 7r^(yl) (^1 + cL2e^2/3n+^ n^(A), observe 
that 

SO that 



1 - no (U„eA{r/. > n+}) " 1 - E^eA^ > ' 
Thanks to monotonicity and Proposition 3.9, for any f G A we have 

(where we took n^/2 to be an integer to simplify the exposition), as required. □ 

Having bounded the probability of exceeding a certain height in the no fioor setting, we can 
now quantify the entropic repulsion effect and derive an estimate on 7r^('^-y /i). 
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Proof of Proposition 3.6. Thanks to Lemma 3.2, it suffices to prove the analogous estimates 
for the measure 11 with no ceihng. 

By Remark 3.8, the conditional distribution of ri\^^ given "^^^/^ is equal to n^^(- | r]\^+ ^ h) 

which stochastically dominates H^^. Hence, 

{v\a, > I '^^,h) > nJ(^(7?rA, > 0) ^ n > 

?;6A^ 

where the second inequality follows by the FKG inequality, the third follows by monotonicity of 
removing the floor and the final inequality by symmetry and Proposition 3.9. Therefore, 

4/3/1 



n^(^rA, > I -r^,/,) > exp (^-2c|A^|e"4^'^j , (3.10) 

since | V (1 — x) ^ exp(— 2x) for x ^ 0. With defined as in (3.6), on the event that 7 is 
an /i-contour and f?(A^) > we have Trj ^ and Il^(Tr]) = e^^^^Il^{r]). It follows from this 
bijection that 

1 ^ ^AiTv) = e^l^'no (r^U^ > , 



■1 

'y,h 



^ exp (^/3|7| - 2c|A^|e-^^'^j K{%,h), (3.11) 

with the second inequality by (3.10). Rearranging this establishes (3.3). To obtain (3.4) note 
first that the proof applies unchanged lihg = h for all s, i.e., when the family of disjoint contours 
is of the form {(7s, h)}s^y , in this case yielding 

n ^7.,/.) ^ exp r ^ (_^|^^| +2c|A,Je-^/^^) ) . 

Now take a general family {(7s, /is)}se.y satisfying the hypothesis of the lemma. We proceed by 
induction over the levels of the contours from top to bottom. If /i+ = max^ hg then conditioning 
on nse,y-hs</i+ '^'ia,hs does not affect the conditional distribution of ??(Use^ ft.s=ft.+ ^t) given 
nse,;^' /is=/i+ '^7s,^+ (^^ explained in Remark 3.8). Moreover, given that flse^^is.'is holds then 
T/j^?7 G ns6-^-/is</i+ ''^7'"'^' ^here T^^ denotes the composition of the T^^'^ for all s such that 
hg = that is, reducing the height of every site in Use^-hs=/i+ This implies that 

n '^7.,'^. I n ^7.,/..) ^ exp r (-/5|7s|+2c|A,Je-4^'^+)]. 

The proof is completed by induction. □ 

Proof of Theorem 3.1, Eq. (3.1). It suffices to prove the corresponding bounds for 11. Set 
h = H — k and Sh{r]) = {v G A : rjv = h}. For each A C Sh^r]) we can define Ua : — )■ given 
by 

{Uav)v = < „ ^ . 

V £ A. 
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To measure the effect of Ua on the Hamiltonian, observe that Ua is equivalent to incrementing 
each height by 1 followed by decreasing the sites in A by /i + 1. As such, this operation increases 
the Hamiltonian by at most \dA\ + 4(/i + 1)\A\ and so altogether 



Hence, 



uKUav) > exp (-4/3L - 4/3(/i + 1)\A\) n^(r?) . 



J2 KiUAV) ^ exp(-4/3L) (l + e''^^'^^'^) H" (r?). 



^ exp (-4/3L + ie-^^('^+i)|5,(r?)|) n^(r?) , 



since e"'^^^'^"'"^^ ^ 1 and {1 + x) ^ e^/^ for ^ x ^ 1. By construction we have Uai] ^ Ua'T] for 
any A ^ A' with A, A' C Sh^rj). In addition, if A C Sh{ri) and A' C Sh ^rj') for some r] ^ rj' then 
C/a?? 7^ Ua'T]' (thanks to the fact that one can recover A from UaV — the sites at level — then 
proceed to recover rj). We can therefore conclude that 

1^ E E 

V-\Sh{v)\ >e--^PkL2 ACSh{v) 



and so, for k ^ 1 

Ki\^h{v)\ > e-'^^^L^) ^ exp (4/3L - \e^P^-^^L 



^ iexp(-e^^L 

where the last inequality holds for any A; ^ 12. A union bound over all k ^ 12 now holds at the 
cost of increasing the pre-factor of 1/2 to 1, as desired. □ 

Proof of Theorem 3.1, Eq. (3.2). As above, we prove the corresponding bounds for H and 
the result for vr will follow from Lemma 3.2. Let /i2 < 2.68 be the connective constant in I? and 
set 



h = H + 



7:5 log 



4/3 Vl-log/i2, 

where Co > is the absolute constant from Proposition 3.6. By the isoperimetric inequality in 
we have |A^| ^ (L/4)|7| for any contour 7 in an L x L box A. Plugging these in (3.3) gives 

Cr^,,,) ^ exp (-/3|7| + Co(L/4)|7|e-^^^') ^ exp(-0|7|) (3.12) 

where 

= /3 - i(l - log ^2) ^ 1 - 3 (1 - log /"2) > log ^i2 

by our hypothesis that /3 ^ 1. 

Now define the random set ^ by 

=2/ = ^(7/) = I7 : 7 is an /i-contour of rj of length I7I ^ log^ L} 

and let be the result of omitting nested contours from ^ : 

s/q = £^o{ri) = £/ e : A^ C A^ for some 7 € 
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For any collection of contours A let also 



Ea= < [j {v e Aj :r]y^ h + k} 
I 76A 



and observe that = E^^ since U{^7 : 7 G =2^^} = U{^7 • 7 ^ M)}- We thus have 

no {E^) = Y,K (Eao I =2^0 = ^0) n^(^o = ^0) 



Ao 



(3.13) 



where 



D ^ h+k • 



iiSA-, 



Conditioned on £/q = Aq, monotonicity enables us to increase the values along A;^ for every 
7 € to /i — 1 while possibly only increasing the probability of the event E^^ , and by doing 
so the variables : 7 G £/q} become mutually independent. 

Fix 7 G If r/i. ^ /i + for some v £ A-y this gives rise to a sequence of nested j-contours 
for j = h + 1, . . . ,h + k surrounding v, and by Proposition 3.6 the probability for a given fixed 
such sequence ipi, ... ,ipk is at most 



exp 



-4/3{/^+i)^ 



) ■ 



However, the fact that Y.j \A^.\e-^^^^+^^ = 0{L^^ log^ L) shows that the area term in this esti- 
mate is negligible, hence the same argument used for proving the upper bound of Proposition 3.9 
(in the no floor setting) yields that, for some absolute c > and every f G A^, 

{ri^,^h + k\ %^h) ^ cexp(-4/3fc) . 

In particular, 
and so 



E 



) / > -^7 
76A0 



nX(-l^7,h)['^7] ^ |A^|cexp(-4/3/£), 



^ cexp(-4^A:) ^ |A^| ^ cexp(-4/3A;)L2 . 
76A0 

The variable "3^ = YlyeAo ^1 therefore a sum of |^o| ^ independent variables, each of 
which respects the bound \S^y\ ^ |A^| ^ log^ L with probability 1. Since /3 ^ 1, for any 
k ^ ilog(4c) we have Ejjo (. ipi^^^^c^^ ^-j [^?^] ^ ^e~'^^^L'^, and applying Hoeffding-Azuma now 
gives 

nO ^ le-^^'^L^ I = ^0) ^ exp (-^e^^^'^L^ log^^ ^) . 
Together with (3.13) we finally get 

H^ {E,^) ^ exp (^-e-4/5'=L2-°(i)) . (3.14) 

Having accounted for this probability in the inequality (3.2), we are left with the problem of 
handling the contribution of long contours, namely those whose length exceeds log^ L. 
Set 

^ = ^{r]) = I7 : 7 is an /i-contour of rj of length I7I > log^ L} . 
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We have shown in (3.12) that, for some 6 ^ Oq with a fixed > log ^2 ^nd any given contour 7, 

Uli'^^^f,) ^ exp(-0|7|). 

By the same argument (appeaUng to Proposition 3.6, this time to the more general bound (3.5)), 
if, for some m = m{L), one considers contours 71, . . . ,7m with disjoint interiors {A^.}^]^ and 
individual lengths all exceeding log^ L, then 



ni(n^7.,/^) < exp(-0f;|7.|) 



■i=l ^ ^ i=l 

By enumerating over the length of each contour 7^, then selecting its origin and a self-avoiding 
path for it (the number of options for the latter being counted by Ri^ii), we see that 

(m \ m 

u u n^^^n E ^^d^-^'^''- 

The relation between 9 and log/U2 suffices to eliminate R\ry-\ while still retaining a factor of 
exp(— |7i|) for some absolute c > 0. The fact that is super-logarithmic now eliminates 

the pre- factor, as well as the additional enumeration over m itself (another polynomial factor). 
Altogether, 



ni( E 1^1 ^ ^ exp (^-ce-^^'L 

for some absolute c > 0, and in particular (via the isoperimetric inequality |A^| ^ (L/4)|7|) 

n^('| U {t; G A^ : r/„ ^ /i} I > \e"'^f^^L^\ ^ exp {-ce-'^^^L^ . (3.15) 

Together with the aforementioned bound on n^(£'_5/), this completes the proof. □ 

4. Lower bounds on equilibration times 

4.1. Proof of Theorem 1: lower bound on the mixing time. 

Set 

h = H-K 

where K is the constant from Theorem 3.1, and define 

e = jr? : #{x E Ai : ^ /i + 1} ^ ^L^ 

Note that, since exp(— 2/3i4r) ^ ^, Eq. (3.1) of Theorem 3.1 implies that 

4(^) = 1-0(1). (4.1) 

Hence, if rg denotes the hitting time of the set it will suffice to show that for a sufficiently 
small constant c > 

minP'?(rB < e"-^) = 0(1). (4.2) 

For this purpose we observe that B is an increasing event so that, 

minP''(rB < e^^) = P^(tb < e"^) ^ P^(tb < e^-^) 

for any initial law v. 
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We now choose v as follows. Take b G (0, ;j) to be a sufficiently small constant so that in 
terms of the constant Cq from (3.3) 

2 

5 < ' 



(/3 - log /^2)^ exp (-4/3(F - /i + 1)) 



where [12 is the connective constant in Tl? . Rearranging the above condition gives 

A := \/^(Co/4) exp (4/?(i/ - /i + 1)) < /3 - log ^2 • (4.3) 

Then we take as starting law v the conditional measure 7r^(- 1 A) where A is the event that there 
exists no /i-contour 7 with area exceeding i.e., 

7:|A^|><5L2 

In the sequel dA will denote the internal boundary of A defined by 

dA := {r\ ^ A: p{r], rj') > for some r]' ^ A} 

where p{-, •) is the transition probability of the dynamics. Let tqa be the hitting time of dA. 

Notice that, up to time tqa, the Glauber dynamics started in ^ \ dA coincides with the 
reflected Glauber dynamics in A whose reversible measure is precisely u = 7r^(- 1 ^4). Therefore 
a simple union bound over times t E [0, e'^^] gives that 

^ e'^ {iy{dA) + iy{B)) . (4.4) 

Define now 

7:|A^|>i5L2 

Notice that dA C A\A since at most four distinct /i-contours can be combined by the modifi- 
cation of a single site. Therefore 

.OA) - ^ < 

We next claim that 

(4 5) 

for some constant ci = ci(/3). Indeed, suppose that 7 is a contour such that |A^|/L^ G {^6,6). 
As in the proof of Proposition 3.6 (see formula (3.11)) and with defined as in (3.6), 

ttHA) ^ ^a(^7^) = e^'"' tA, > I A, ^,,h)A{A n %^h) 

'?6^,r;rAT,>0 

where we used the fact that T^ry £ A if r/ G ^ n "^^.h,. Next we observe that, thanks to (3.10) 
(which holds with identical proof also for vr^^), 

^(r^rA, > I ^,-^7,0 = > I '^-y,/^) ^ exp(-2c|A^|e-^^'^) 

to yield 

4 {"^^^h I A) ^ exp (-/3|7| + ColA^I exp(-4/5/i)) . (4.6) 
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The isoperimetric inequahty in gives that |A^| ^ I7P/I6 for any 7, so that, by the above 
choice of parameters, any contour 7 with area less than bl? satisfies 

C^\k,\e-'^^ ^ Co (Vs^yWN) (^J«"'^' 

< ^/5(Co/4)|7|e^^(^-^+i) = A|7| 

where A is given by (4.3). Hence the r.h.s. of (4.6) is smaller than e"^'^"^^''^' . A union bound 
over 7's with |A^| > ((5/5)L^ then proves (4.5). 

In conclusion the first term in the r.h.s. of (4.4) is o(l) if c < ci. We now examine the second 
term ^{13) and we proceed as in the proof of Theorem 3.1. First, we claim that for any short 
/i-contour 7 and f E A^, where "short" means of length smaller than log^(L), we have 

(??. ^h + i\ ^^^h) ^ \ , (4.7) 

Indeed, if = {7' : v E A^/ C A^} then an application of (3.5) from Lemma 3.7 shows that 
vri (r?„ ^ + 1 I %,h) ^ 5^ ^A (^7',/^+i I ^7,/^) ^ ^ I 

^'^si 7' Si/ 

for /? large since, as usual, the number of contours 7' E iz/ of length k is at most kyl^ (using the 
fact that each of these crosses the horizontal line to the right of v within distance at most k). 
To transfer this estimate to the setting of a floor, observe that by Remark 3.8, 



(7?, ^ /i + 1 I ^^^h) = — ^ . (4. 



We have just established that the numerator is at most 1/8, whereas by monotonicity the 
denominator is at least 

< (^fA, 0) = 4, (^rA, ^ - /i) ^ 1 - ce-'^^^+^M,\ 

thanks to Proposition 3.9 (with the same constant c > appearing there) and a union bound 
over the sites of A^. The fact that |A^| ^ |7|^ = 0(log^ L) shows this last term is 1 — L~'^'^°^^\ 
hence the effect of the denominator in (4.8) can easily be countered by a factor of 2, thus 
establishing (4.7). 

With inequality (4.7) available to us, the very same concentration argument leading to (3.14) 
applies again here to imply that 

"a (E/i^ E A, : r?,, ^ /i + 1} ^ ^ e~^-^'^"''' (4.9) 

for some constant C2 > 0, where the summation is over every short /i-contour 7. Similarly, 
following the same steps leading to (3.15), we get that 

K (E/{^ e A^ : ri^^h + l}-^ ^ e-'^^^ (4.10) 



for a suitable C3 > 0, where sums over every long /i-contour 7, i.e. such that I7I ^ (logL) 



2 



(in this case, the analog of (3.12) for /i-contours of area smaller than SL"^ holds if 6 chosen 
small). Finally, Lemma 3.2 translates the statements on 11^ into the analogous bounds for vr^. 
In conclusion the second term in the r.h.s. of (4.4) is o(l) if c < min(c2,C3), as required. □ 
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4.2. Proof of Theorem 2: lower bound on Tq. Here we prove that P'^(ra ^ e^^") — 1 as 
L — >■ oo where, we recall, 

Qa = {r] such that #{x £ Al ■ > aH{L)} > ^|Al|} 

and Ta is the hitting time of ^la- We proceed as in the proof of (4.2) but now the height h is 
chosen equal to aH{L) — 1, so that e~^^^ ^ exp(8/3)L~°, and the set A is defined by 

7: |A^|>(5L2" 

Here 5 is a small constant such that, for |A^| ^ 5L^": 

ColA^le"^^'^ ^ Co (\/^vW7l6) e"^'^'* ^ AI7I 
where A is analogous to (4.3). As in (4.4) we get 

F'^iTa < 6^=^") ^ F'^iraA < e^^") + F^r^ < e^^" ^ T9a) 

^ e-^^" (i/(M) + 1/(0,)) . (4.11) 

Exactly the same arguments behind (4.5), (4.9) and (4.10) now show that the r.h.s. of (4.11) is 

0(1). □ 



5. A BOUND USING PATHS AND FLOWS 

5.1. Proof of Proposition 2.3. Let A := {1, . . . , L} x {1, . . . ,m} and Q, := f^An+- We intro- 
duce the canonical paths 7(??, ?/) from ij to rj' for every rj,?]' £ Q. Define the diagonal lines in 
Al = {1,...,L}^ 

Ri = {x £ Al ■■ X2 = xi + L - i}, i = l,...,2L-l (5.1) 

and let TZ denote the collection of the Ri. Number the sites in A following the lines Ri, . . . , R2L-1, 
so that each line is read from south-west to north-east; at each site x move straight from r]^ 
to r]'^ by taking \r]x — unit steps. Note that since all heights satisfy ^ r/a; ^ one has 
I7I ^ lAIn"*". If e = ((7,(7^*'^) is an edge of a path, with x^ G , define A as the set of 2; G A 
such that X < x^ and B the set of x > x* (w.r.t. to the order introduced above). Here o"^*'^ 
denotes the configuration which coincides with a except that the height at x* is changed by ±1. 
Then by direct inspection one finds that for any r/,r/' G such that 7(f?, ??') 3 

7r(77)^(7?') ^vr(a)7r(a*)exp(6/3 ^ |r?x-r?^|) (5.2) 

xei?i, nA 

where o satisfies (ta = Va^ '^B = Vb, while a* is the configuration obtained by setting (t\ = ija, 
(^*B = ^'b- Here and o"*^ = ± 1 are assigned according to the choice of e. The crucial 
observation is that, given e, the map from (?/, ?/) (such that e G 7(7?, 77')) to (c, a*) is an injective 
one. In particular, this implies: 

— r-T |7(^,^')k(r?)vr(?7')lee7 ^ l^l'^'^exp(6/3n+m). (5.3) 

Note also that the inverse of the smallest non-zero one-step transition probability for our chain 
is |A| exp(4/3n^). We apply then (2.10) to obtain that the inverse spectral gap of the SOS 
dynamics is upper bounded by clApn"*" exp (7/5n'^m) and (2.11) follows. □ 
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5.2. Proof of Theorem 2.4. For every G il, let 71 be a path of length T starting at ^ 
and let 72 be a path of length T starting at Write r/, r/' for the corresponding endpoints. Let 
7c be a path from to r]' (to be specified below) which depends only on r/,r/' and not on 71,72- 
Call 7 the concatenation of 71, 7c, 72) where 72 is the path 72, inverted in time. Note that 7 
connects to If rj,?]' £ G then we let 7c be the path 7(7/, r/') which appears in the statement 
of the theorem (recall that it stays in the set G) and define 

7r(g)F^(7i) vr(Opg'(72) 
P«(X(T) G G) P«'(X(r) G G) 

Otherwise, set 0(7) = and we do not need to specify 7c in this case. Here P^(7i) is the 
probability that the process {X{t))t started at follows exactly 71 up to time T, and similar for 
P^ (72). Note that for fixed G Q, X]^.^^^/ 0(7) = 7r(^)7r(^') where the sum is over 71, 72 
for fixed ^, 

Therefore, viewing the path 7 as a collection of oriented edges e = {cr,a') and letting Vg/ = 
f{a) — f{(y')-, we have 

Var(/) = ij;^(07r(O(/(6-/(e'))' = ^E E <i)[Y.'^ ef)' (5.4) 

5,e 5,5' 7:5-5' £67 

^^E E a(7)(A(/) + ^7(/)), (5-5) 

5,5' 7:5-5' 

where 

A(/) = i7iiE(^£/)' + i^2iE(^£/)' 

6671 6672 

B,{f) = \^\Yjyef? 

667 

and in the inequality we used Cauchy-Schwarz. Now we use the fact that the Dirichlet form 
which appears in the definition (2.5) of the spectral gap can be written as 

f(/):=4(/(/-P)/) = i E ^(^M^,^')(Ve/)^ 

e=(cr,o-') 

Recall that Pmin denotes the smallest non-zero one-step transition probability, and observe that 

5,5' 7:5~5' 5 71 V V 7 ) gg^j 

^^^—£U) sup (^(a)-i)EE^(^)^^(^i)W' 

where we used the fact that |7i| = T. 

Let P denote the law of the stationary process (started at equilibrium vr). From a union 
bound, one has 

EE^(^)'P^(^i)W = G [0,r] : X{t) = a, X{t+l) = a') ^ T7T{a). 

5 71 

It then follows that 

E E a(7)^(/)^ £:(/). 

5,5 7:5~5 
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As for the second term in (5.5), using stationarity of vr one has that the sum of ■7r(^)P^(7i) over 
all ^ and paths 71 of length T which connect ^ to r] gives Tr{r]), so that (with the definition 
(2.12)) 

E E oh)B,(/)<|i (v./)^.(.)p(.,,o E 

^^W{G)£{f). 

Going back to (5.5) and to the definition of spectral gap one immediately gets (2.13). □ 

6. Upper bounds on equilibration times 

6.1. Proof of Theorem 2: upper bound on Tq assuming Theorem 1. Here we prove 
that P'^(ra ^ e*^^") — )• 1 as L — )• oo assuming T^ix =^ e*^^. The latter estimate will be proven 
afterwards. Let us partition the box into non-overlapping squares Qi of side CL" with 
C = exp(4/3K) where K is the constant appearing in Theorem 3.1. By monotonicity the Glauber 
dynamics is higher than the auxiliary dynamics in which each square Qi evolves independently 
from the others with boundary conditions on dQi. Using the assumption Tmix ^ e'^^ and 
independence, it is standard to check that the mixing time of this auxiliary dynamics is not larger 
than e^'^^" and therefore, at time T = e^'^^" , all the squares Qi are close to their equilibrium (in 
total variation) with an exponentially small error. Theorem 3.1 implies that in each of them the 
density of vertices higher than 

H{CL'') -K = aH{L) 

is larger than 1 — e(/3) with probability exponentially close to one. In conclusion, apart from 
an exponentially small error, f^{Ta > e^'^^'^) is bounded by the probability that for some i the 
square Qi has a density less that 1 — e{(3) of vertices higher than H{CL"') — K. Thus, a union 
bound suffices to conclude the proof. □ 

6.2. Proof of Tmix < for n+ = log L. To prove the upper bound on T^ix in Theorem 1, 
the crucial point is to give the proof for n+ = log L, so we assume this is the case in this section. 
The general case logL ^ n+ ^ L can be then deduced via very soft arguments; see Section 6.3 
below. 

For reasons that will be clear later, first of all we modify the SOS model by considering the 
Boltzmann factor exp[— /3?^^^ + /] instead of exp[— where / is the external field term 

^ n+-H n+-H 

f=Y^^fy^ with fy= ^ fy,j-= Yl ^j'^Vy^H+j, (6.1) 

s/eAi j=i j=i 

with H = H[L) defined in (1.3) and Cj = exp(— /3j). One changes the partition function 
accordingly. We call vr^'^ the corresponding equilibrium measure with ceiling at n"*" = log L and 
floor at 0. Moreover, we will consider the Glauber (heat bath) dynamics associated to tt^'^. 

Remark 6.1. Note that, if the b.c. are zero then the extra term f in (6.1) will not drastically 
change the global equilibrium properties, since it tends to depress the heights that exceed the level 
H (and having rjx ^ H + 1 is already an unlikely event, for (3 large). More precisely, f equals 
the constant (|A/,|/L)^jCj plus a (negative) random term which one could prove, by refining 
the estimates of Section 3, to be of order L x exp(— c/3) for a typical configuration (and therefore 
not extensive in the area of Al). 
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The reason for modifying the equihbrium measure in such a pecuhar way is explained after 
Theorem 6.12. 

Lemma 6.2. The ratio A of the mixing time of the original system over the mixing time of the 
system modified as in (6.1) satisfies for L large 

^ A ^ e^. (6.2) 

Proof. Going back to the definition (2.5) of the spectral gap, it is easy to see that the ratio A 
of relaxation times satisfies 

e-4|/U ^ A ^ e^\f\^ 

with / as in (6.1); see e.g. [39, Lemma 13.22] for such standard comparison bounds. Note that 
|/|oo = 0(Le-^) if/3is large enough. Then (6.2) follows from the comparison (2.7). □ 

Therefore, it is enough to prove Theorem 1 for this modified model. We denote its mixing 
time as Tyii^{L). It is important to realize that the Glauber dynamics for this modified SOS 
model is still monotone (in the sense of Section 2.2) and that the FKG inequalities are still 
valid. This is because / is the sum of functions of a single height r]x- Therefore we can ap- 
ply all the monotonicity arguments we need (including the Peres- Winkler censoring inequality. 
Theorem 2.2). 

Definition 6.3. For A; G N and o, A > 0, we define the inductive statement J-^ := J-k.a,A ■ for 
every L the mixing time satisfies 

where log^'''^(2;) := max (1, log (log ...(x))) and log (log ...(x)) is the logarithm iterated k times. 

Theorem 6.4. Fix /? ^ /3o for some large enough constant and = log L. Then =^ J'k+i 
provided that a = 4 and A = for some sufficiently large C . 

Proof of Theorem 1 given Theorem 6.4- For k = 1 the statement J-i follows at once from the 
"canonical paths argument", Proposition 2.3 (with o = 3 and A = bp, b some explicit constant). 
Notice that Proposition 2.3 applies with no change to the modified model with the external 
field. Then, apply the theorem until log^'^-'(L) = 1. At that point we get the desired exponential 
mixing time upper bound. □ 

For the proof of Theorem 6.4, we need some notation. Recall the definition (5.1) of the 
diagonal lines Ri. Define C J^l as the set of configurations r] such that, for every R ^TZ, 

(with [x]"*" = max(x,0)), C JIl as the set of configurations rj such that, for every R £ TZ, 

Y,[H - r,,]^ ^ LI 

and finally Gf C 0,^ as the set of configurations r] such that, for every R £ TZ, 

Y,\H -rjx\i^LL (6.3) 
xeR 

Let also 

£(A;,L) :=SlogW(L) + ^log^ (6.4) 
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with B a constant to be chosen sufficiently large (independently of /3) later, see discussion after 
(6.7) and (6.20). 

Lemma 6.5. Assume with a = 4 and A = 40BP and take Ti > e^^. Then 

P(r?^(Ti)GG-,_^,^^))^f. 

Lemma 6.6. Assume J-k with a = 4 and A = 40Bf3 and take T2 > e^^^ with the same B as in 
(6.4). Then 

p(^^(r2)GG+^^^p^|. 

Note that ^(00, L) = B + 1/(4/3) log A is just a large constant. We will actually see that, in 
both lemmas, the constant | can be replaced by 1 — o(l) where o(l) vanishes for L — t- 00. We 
refer to Section 6.4 and Section 6.5 below for the proof of Lemma 6.5 and Lemma 6.6 respectively. 

Proof of Theorem 6.4 given Lemmas 6.5 and 6.6. Thanks to monotonicity, to Lemmas 6.5-6.6 
and the fact that ^ C G^^^-^ we can set T^" := max (e^'^,e^^'^) and obtain that 

mmP(77«(r-ll) G ^2,(^+1,^)) ^ \. (6.5) 

This is based on the fact that, if rj^ ^ 7] ^ rf' and E ,rj^ G then rj G G^+^z. Just write 

|r?. -H\ = [r?, - H]+ + [H - r?,]+ ^ [rg - H]+ + [H - r?^] + . 
At this point we need the following consequence of Theorem 2.4. 

Proposition 6.7. Let a,£ > and T be such that ¥{r]'^{T) G Ge) ^ a for all initial configura- 
tions C- Then there exists a constant c = c(a,/3) such that 

T^el{L) ^ c [exp (15/3£L) + L^^ T^l . (6.6) 

Proof. Theorem 2.4 gives 

J-REL ^ ~ 1 



^ \Pmin OL 

where in the definition of W{G() we choose the canonical paths introduced in Section 5.1. We 
know that the invcrsG of the niinimajl tra,nsition probcibility Pxnin is of order |Ai| exp(4/3n+). 
Also, from the proof of Proposition 2.3 and the definition (6.3) of G^, we see easily that 
W{Gi) < exp(15/3^L) and then the claim fohows. □ 

Proposition 6.7 (applied with T = T'^", £ replaced by 2£(/c + 1,L) and recalling that n+ = 
logL), together with (6.5) and (2.7), implies that 

Tmi.{L) ^ c'iP)L' ( L^^iT^y + e30i?/3ilog(^+i)(L)+8LlogA\ _ ^g_7) 



If one chooses A = 40BP (and B large but independent of k, (5) then the r.h.s. of (6.7) is smaller 
than L4exp(ALlog(''+^)(L)) for every L and the claim follows. □ 
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6.3. Proof of Tmix ^ e.^^ for logL ^ n+ ^ L. Once we have the statement for n+ = logL, 
proving it for logL ^ n"*" ^ L is quite easy, so we only sketch the main steps. Thanks to (2.8), 
it is enough to prove that 

- 7r|| ^ (6.8) 
ll/i^ -7r|| ^ (6.9) 

for some t = exp(0(/3L)). Here we write vr instead of vr^ for simplicity. We first note that, if 
7r,7r are the equilibria with ceiling at > logL and at logL respectively, then 

||7r-7f|| ^ L^'o^'^) (6.10) 

with co(/3) that diverges as /3 — )• oo. Indeed, to feel the ceiling there must be some x such that 
77a; ^ logL and this has probability at most c\Al\ exp(— 2/3 log L). This can be seen as follows. 
By monotonicity lift the b.c. from to (logL)/2. In this situation, the probability that the SOS 
interface reaches either height or logL is 0{\AL\e~'^^^°^^), as follows from Proposition 3.9 and 
a union bound, cf. the proof of Lemma 3.2. 

As for (6.8), from Theorem 2.2 (applied with k = 1, ti = t,Vi = Ai,ai = 0, 6i = logL) 
we have Wn^ — n\\ ^ — 7r||, with jit the law of the evolution fj{t) with ceiling at logL. 
Since we proved in Section 6.2 that the mixing time of the dynamics fj(t) is exp(0(/3L)), if 
t = exp(c/3L) with c large one gets from (2.4) that — 7r|| = o(L~^) and therefore — 7r|| = 
0(L-co(/3)) + ii^u _ ~ II ^ o(L-4) if /3 is large enough. 

As for (6.9), assume for definiteness that n"*" is a multiple of logL and let 

i — 1 2n+ 

hi = n~^ —logL, i = l,...,M := 

2 log L 

Let us apply Theorem 2.2 with k = M, Vi = Al, ti = iexp(c/3L) with c large enough, 5j = hi 
and ai = hi — logL. Let us also call Ui the event that ^ r/^, ^ Oj + ^ logL for all x ^ Ai. 
Note that, for the associated modified dynamics ??(t), in the time interval < t ^ ti = exp(c/3L) 
the floor is at height ai = n"*" — logL and the ceiling at height bi = n"*". Therefore, if c is 
chosen large enough, at time exp(c/3L) the system is within variation distance say from 
the equilibrium with such floor /ceiling and in particular, except with probability smaller than 
^-co(/3)^ the conflguration is in Ui (the proof of this is very similar to the proof of (6.10) above). 
If r/(ti) G Ui, then in the second time-lag {ti + 1, . . . ,t2} the situation is similar, except that 
the floor is now at 02 and the ceiling is at 62 (note that if instead fj{ti) ^ Ui then some heights 
are frozen forever to values larger than oi + ^logL and the dynamics ?7(i) will not approach 
equilibrium). The argument is repeated M times with the result that (via a union bound on i), at 
time tM = exp(0(/3L)), the variation distance from equilibrium is smaller than ML"'^"^^) ^ L~^ 
and the proof is concluded. □ 

6.4. Rising from the floor: Proof of Lemma 6.5. We will make a union bound on Ri G TZ, 
i.e. on i = 1, . . . , 2L — 1. We want to upper bound 

^[F-r?^(ri)]+^Lf(fc + l,L)) . (6.11) 

Reader's Guide 6.8. 

In principle the argument is very simple. Around every point x € Ri one would 
like to consider a square Qx of side L^ := L/(Alog^^^(L)). By monotonicity the 
quantity [H — t]);! (Ti)]'^ appearing in (6.15) gets larger if we fix to the heights on 
dQx- Prom the assumption J^k we know that the mixing time in Qx, with zero b.c. 
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on dQr^, is of order exp(^Lfc \og^^\Lk)) exp(L) ^ Ti exp(2L). Thus, at time 
Ti the dynamics in Qx is essentially at equilibrium (w.r.t. zero b.c. on dQx), so 
that r]x ~ 1/(4/3) log Lfc ^ H - (1/4;S) log(''+^)(L) w.h.p. By taking the constant 
B appearing in (6.4) large enough we can make i{k + 1, L) ^ (1/4/3) log''^^"'^^(L). 
As a consequence, the event in (6.11) describes a very unlikely deviation. 

In practice, the proof is considerably more involved, in particular because the 
size of the squares Qx has to be chosen as a function of x (cf. (6.16)) in order to 
guarantee that Qx is fully contained in the original domain A^. 

Set ax '■= logd(x), where d{x) is the distance of x from the boundary of A^. One has 

[H - Vx{T,)]+ ^ [ax - r?^(ri)]+ + - H\. (6.12) 
Now, there exists Ci such that for every Ri G TZ one has 

Y^\ax-H\i^^L. (6.13) 

By the way, this is the reason why we defined the lines Ri as in (5.1): if Ri were parallel to 
the coordinate axes and too close to the boundary of A^,, then (6.13) would be false. To prove 
(6.13), suppose without loss of generality that the diagonal line under consideration is Ri with 
i ^ L, so that \Ri\ = i. One has 

E \ax-H\=iH- ^^i(logL-i J]logd(x)). (6.14) 

xGRi xdRi xdRi 

If i is even 

i/2 

log d{x) = 2 E log k = i log i + 0{i) 

x^Ri k=l 

and a similar argument takes care of the case where i is odd. Therefore, 

C"L 



V \ax - H\ ^ — L 

I ^ ' ^ 4/3 
xeR, ^ 



'-\og{L/t) + C'^ 



/3 



for some constants C", C" > independent of i, f3,L. 

Let us go back to estimating (6.11). It is clear that the x such that d{x) ^ L/logL can 
give altogether a contribution to "^xi^x — "^x (^i)]"*^ which is at most 0{L). Then, let Ri be the 
subset of Ri such that d{x) > L/logL. We can conclude that it is enough to estimate 

J2[<^x-VxiTir>mk + l,L)-C"') 

\xeRi 

^F(^[ax- VxiTi)]^ > ^^(^ + hL)]. (6.15) 
\xeRi I 

Now for every x ^ Ri define a (diagonal) interval Ix C Ri, centered at x and of length 

\Ix\ = ^min ld{x), ^ ) . (6.16) 

Note that the minimal | Ix \ is of order L / log L and the maximal one is at most L/(2Alog('')(L)) (it 
can be much shorter if \Ri\ <^ L/log^^\L)). Note that condition (6.16) guarantees that around 
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each Ix one can place a square Qx of side rux = '2\Ix\ and fully contained in A/,. Considering all 
the possible i ^ L and the different intervals Ix,x 6 Ri, the number of such intervals is trivially 
smaller than |Ax,|. Therefore, observing that Ri can be covered by (possibly overlapping) such 
intervals Ix of total length at most (3/2)|i?j|, it is enough to prove 

rlY.^ay-v'^in)]+;,\^i{k + l,L)] (6.17) 

for every such interval and then apply a union bound to get that the r.h.s. of (6.15) is o(l/L), 
so that after summing over the index of Ri the probability in (6.11) is still o(l). 
It is easy but crucial to check that 

T^^cTy-v^inr < Y.^H{mx)-v^{T,r + ^\Ix\{logA + log('+'\L)) (6.18) 
y&Ix y&Ix 

where of course, as in (1.3), H{mx) = 1/(4/3) log mx is just the typical equilibrium height of the 
SOS interface in the center of the square Qx with zero boundary conditions on dQx (here, for 
lightness of notation, we forget the integer part in the definition of H). Indeed, since nix ^ d{x) /2 
one has for y & Ix 

\ay - H{mx)\ = ^(logd(y) - log(m^)). (6.19) 

If min((i(x), L/Alog^'^^(L)) = d{x) then the r.h.s. of (6.19) is upper bounded by a constant. In 
the opposite case, it is bounded by 



1 



log L — loe 



L 



and (6.18) follows. Therefore, it is enough to bound 



^^(log^ + log('^+i)(L)) 




[H{mx)-Vy{Ti)]+ ;^Co\Ix\j ^L-' (6.20) 

for all such intervals, for some Cq independent of /3. We can assume that Co is large (just choose 
B large in (6.4)). 

Monotonicity implies that if we let evolve only the heights inside Qx with 0-b.c. on dQx, 
then the random configuration obtained at time Ti is stochastically lower than the configuration 
obtained via the true evolution (where all the heights are updated). Again by monotonicity 
(the event in (6.20) being decreasing) we can lower the ceiling in the box Qx from height 
n"*" = logL to height logm^ and also replace the pre-factor (l/L) with {l/nix) in front of 
the fields fy,y £ Qx in (6.1): the dynamics thus obtained (that we simply call "the auxiliary 
dynamics") gets stochastically lower. The reason is that the fields fy are decreasing functions 
of r], which tend to "push down" the interface, and l/rux > l/L, so that exp((l/L — l/mx)fy) 
is an increasing function. 

Since we are assuming that J-k holds (with a = 4), the mixing time of the auxiliary dynamics 
in Qx (with 0-b.c. on dQx) is at most 

mxex.p{Amx\og^''\mx)) ^ L^exp(L). 

As a consequence, using (2.4), at time Ti = e"^^ the law of the auxiliary dynamics is within 
variation distance exp(— e^/^) from its invariant measure, call it ttq^, which is nothing but a 
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'Of -Of 

space translation of t^/^^^^ , where we recall that, for a generic L, ■k^-'^ is the equilibrium measure 
in Kl with the field /, the floor /ceiling constraints ^ ^ log L and b.c. at zero. For simplicity, 
for the rest of this subsection, we shift the square Kl so that its center coincides with the origin 
of I?. 

In conclusion, 

and / is a diagonal segment of cardinality |/| = = f^'^xl'^i centered at the origin of I? . Thus 
we need the following equilibrium estimate: 

Lemma 6.9. For any m, if I is a diagonal segment of length \I\ = m/2 centered at the origin 
of , then: 

-KlliB) := vr^'f {jl\H{m) - r?,]+ ^ Co|/| j ^ cexp(-/3m/c) (6.21) 

where c> is a constant and Am denotes the side-m square centered at the origin. 

This will then be applied with m ranging from order L/logL to order L/log^*'^(L) so in all 
cases the r.h.s. is much smaller than and, putting everything together, the inequality (6.20) 
and therefore the claim of Lemma 6.5 follows. □ 

Proof of Lemma 6.9. Suppose this is true for the model without the field /, i.e. for the standard 
SOS measure vr^^^^ of (2.2). Then, the same estimate follows (for /3 large, with c replaced by 

c/2) for vr^'-^. This is so because, uniformly. 



m 

for some c' independent of /3. To show that 7r^^(;B) is small, one first proves that 

TT^'") {B) ^ exp(-(Co/4)/3m) (6.22) 

say for every |/| of size between and |m, where we recall from Section 2.1 that vr^^™^ is the 
SOS measure without floor/ceiling and boundary conditions at height H{m). This is based on 
Peierls-type arguments and the proof is relegated to Appendix D. 

We conclude the proof of Lemma 6.9 assuming (6.22). Define Aj, i = 1, . . . , m/2 — 1 to be the 
boundary of the square of side m — 2i centered at zero. Let Ei be the event 

Ei = [Y [H{m) - r?J+ ^ 5m} (6.23) 

for some 6 to be chosen small later. Suppose that at least one of the Ei,i ^ m/10 is not realized, 
and let j be the smallest such i. In that case, we look at the vr^^ -probability of B, conditionally 
on the configuration of r] on Aj. For all x € Aj, if ijx > H{m) we can lower it to H{m) by 
monotonicity (the event B is decreasing). If instead r]^ < H{m), we still change jy^, by brute 
force to H{m): the price to pay is that in the final estimate we get a multiplicative error 

exp lc(3Y[H{m)-rj,]+] ^ e'^'^ 



30 



P. CAPUTO, E. LUBETZKY, F. MARTINELLI, A. SLY, AND F.L. TONINELLI 



for some explicit c (independent of (3 and 5). What we get is that, conditionally on j ^ m/10 

i^j is not realized, the vrj^^ 



being the smallest index such that Ej is not realized, the vr? -probability of B is upper bounded 



by 

gC/35m^HM^ (^|o ^ ^ ^ log m) ^ e'^^'^T^f^X, (B\r^ ^ log m) (6.24) 

where the inequality is just monotonicity. Notice that vr^*-'"^ (r? ^ log m) is large (say, larger 
than 1/2, cf. Proposition 3.9). Then, we can apply (6.22), since the interval / we are looking at 
is of length m/2, so that certainly \{m — 2j) ^ \I\ ^ |(m — 2j) and we get that the r.h.s. of 
(6.24) is upper bounded by 

exp(c/3<5m - (Co/4)/3(m - 2j)). 

At this point it is enough to choose 6 small enough, for instance S = Co/ (20c), to conclude 
(recall that j ^ m/10). 

Next, we have to show that 

Kjr^TJl'EA (6.25) 



is very small. Indeed, that event implies that 



^ [H{m) - 7?^]+ ^ JmVlO = ComV(200c). (6.26) 
Write 

^ [H{m) - = J2 (6.27) 

where Ni^. is the number of points where [H{m)—rix]~^ = k. From Theorem 3.1 we know that there 
exists some integer K such that Nj^ ^ m'^e~^^^, except with probability exp(— ?nexp(/3A;)), for 
k ^ K. Then, except with probability of order exp(— c/3m) one has ^ ^ kN^ < Com^/ (200c) 
if Co is chosen large enough (recall that, as discussed after (6.20), we can assume that Co is 
large) . □ 

6.5. Falling down from the ceiling: Proof of Lemma 6.6. This is the part which requires 
the more subtle equilibrium estimates. Let T2 = exp(c/3L) where c will be determined along the 
proof. We want to prove that 

F(r?^(r2)GG+^^^))>f. (6.28) 

We recall that i{oo,L) = B + 1/(4/3) log ^ is a constant that we can assume to be large. For 
simplicity, we write i instead of ^(00, L). 

Reader's Guide 6.10. 

Ideally the proof would work as follows. At equilibrium, the event Gf has prob- 
ability almost 1, see Lemma 6.11 below (since is decreasing, in Lemma 6.11 
we lift the boundary conditions on OAl from to H' = H + 1, the reason for the 
"+1" being that, in this way, for /? large the floor has little influence on the inter- 
face at the typical height H' .) It is therefore sufficient to prove that at time T2 
the dynamics (with b.c. 0) is close to equilibrium. For this purpose we will apply 
Theorem 2.2 (which is allowed since we start from the maximal configuration □) 
with the following censoring schedule. 

Cover Al with overlapping, parallel rectangles Vi, i ^ M = O(logL), or- 
dered from left to right, with longer vertical side L and shorter horizontal side 
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(L/ (log L)) and such that ViDVi+i is a rectangle L x (L/(21og L)). Now consider 
the "bricks" Bi which have base Vi and height n+ = logL. 

We first let Bi evolve for a time ti = exp((c/2)/3L). This is the SOS dynamics 
with b.c. on the left, top and bottom boundary of Vi, and with b.c. on the 
right boundary. As we justify below, we can pretend that at time ti, the system 
in Bi has reached its own equilibrium. This equilibrium, restricted say to the 
left half of Bi, should be extremely close to the true equilibrium in Al with 
b.c. This can be justified as follows. The b.c. around Vi impose the presence of 
open contours at heights 1, . . . , n"*", with endpoints at the endpoints of the right 
side of Vi. These contours behave roughly like random walks and will stay within 
distance say from the right side of Vi and only with tiny probability will 

intersect the left half of Vi . 

Next, we let B2 evolve for the same amount of time ti, after which a similar 
argument shows that the "true equilibrium" is reached in the left half of V2, 
i.e. on the right half of Vi, and so on. When the M*'* block has been updated, 
the system should be very close to equilibrium everywhere. In practice, there 
are two major obstructions that prevent this strategy from being implemented 
directly and which cause much technical pain. The first has to do with the 
presence of the floor constraint at zero and will be discussed in greater detail in 
the Reader's Guide 6.13 below. The second difficulty can be understood in the 
following simplified situation. 

Take the SOS in a L x m rectangle R, with \/Z <^ m <^ L (for us, R would 
be Vi so that m = L/logL) with b.c. 1 on one of the size-L sides and b.c. 
everywhere else, without any floor/ceiling. There is an open 1-contour joining 
the endpoints of the side with 1 b.c. The probability of such contour 7 can be 
shown, via cluster expansion, to be proportional to exp(— /3|7| + ^'^(7)) where 
the "decoration" term ^ r{'^) is of order I7I times a constant which is small with 
(3. In absence of decorations, 7 would behave as a random walk and it would be 
very unlikely that it reaches distance ^ \/Z from the side with 1-b.c. In presence 
of the decorations, this might in principle fail. Indeed, the decorations depend 
also on how close the contour is to the boundary of R (see Appendix A), and this 
could induce a pinning effect of the contour on the size-L side with 0-b.c. The 
way out we found to exclude this scenario is a series of monotonicity arguments 
which in practice boil down to transforming R into a rectangle with both sides of 
order L. In this situation, since the side with 1-b.c. is very far from the opposite 
side, the "pinning effect" can be shown not to occur. 

To prove (6.28) we couple rp{T2) with a suitable equilibrium distribution as follows. Let A 
be the 2L x L rectangle obtained by attaching a square of side L to the left of the original 
square A^. Let vr^ '■' denote the SOS equilibrium distribution in A with boundary conditions 
H' := H + 1. Such equilibrium measure contains the field /, cf. (6.1) (where the sum now is 
over y £ A and the pre-factor is still 1/L) and floor/ceiling constraints ^ 77 ^ n'^ . One has 

Lemma 6.11. If i is large enough, then 

lim vrf '-^(77 G G+) = 1. (6.29) 

L—^oo 

The proof is deferred to Appendix D. 

Therefore, using that the event is decreasing, (6.28) follows if we prove that there exists a 
coupling of {rj, rp (T2 ) ) , where r] is the restriction to Aj;, of the configuration distributed according 
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to TT^ '-^ , such that 



P(??^(72) ^ 7?) = l + o(l). (6.30) 



To this end, we will apply Theorem 2.2 with exactly the censoring described above. We first let 
evolve the system in Bi for a time-lag ti, with n"*" b.c. on the right side of Vi and b.c. elsewhere. 
Then we let evolve the system in B2, for another time-lag ti. For B2 we have the maximal b.c. 
n'^ on the right boundary, zero b.c. on top and bottom and the b.c. on the left boundary is 
given by the configuration, say ri, inherited from the previous evolution on Bi. We repeat this 
procedure for the other bricks Bi, i < M, with maximal b.c. on the right boundary, zero b.c. on 
top and bottom and the b.c. rj_i on the left boundary; the final brick Bjyj, unlike the previous 
ones, has a zero b.c. on the right boundary as well as on the top and bottom boundaries, and 
b.c. tm-1 on the left boundary. 

We let rj denote the configuration at the end of the above described procedure. Note that 
altogether the time spent is Mti ^ T2 = exp(c/3L). Theorem 2.2 implies that we can couple 
rj and r]'^{T2) in such a way that P(7y'^(T2) ^ rf) = 1. Thus, it remains to prove that (6.30) is 
satisfied with replacing if^{T2). 

The mixing time of a brick is bounded above by exp((c/4)/3L), for a suitable choice of c > 0, 
see Proposition 2.3. Therefore, after time ti the chain is extremely close to its equilibrium in Bi 
with the given boundary conditions. Up to a global error term of order we can thus assume 
that after each updating of a brick, the corresponding random variable is given exactly by the 
equilibrium distribution on that brick with the prescribed boundary conditions (see Eq. (2.4)). 
Let rji denote the configuration after the updating of brick Bi, restricted to the left half of the 
brick, i.e. the brick with basis ¥[ := Vir\ (V^+i)'^. Thus, using monotonicity, it is sufficient to 
exhibit a coupling such that 

P(^i ^ r/i, i = 1, . . . , M - 1) = 1 + 0(1), (6.31) 

H' f 

where rji denotes the configuration t] with distribution vr^ '■' , restricted to . 

To prove the latter estimate we proceed as follows. Let V* denote the portion of covered 
by rectangles Vi, . . . ,Vi, and set V'^ := V{. For i = 0, . . . , M, call A* the rectangle obtained 
by attaching a square of side L to the left of V* (this corresponds to the "rectangle enlarging 
procedure" outlined above), and let ^ denote the b.c. equal to: 

^ I if X belongs to the right boundary of A* 

I H' otherwise. 

Since C A'~^, by monotonicity and a repeated application of the DLR property for the 

measure vr^ , we see that the desired claim (6.31) is a consequence of the next equilibrium 
result. 

Theorem 6.12. For every C > 0, there exists /3q such that for all /? ^ /3q, for all i = 1, . . . , M , 

where \\ ■ denotes total variation of the marginal on A*. 

Reader's Guide 6.13. 

We now explain why (6.33) should be true and why we crucially need the field 
/, which is absent in the standard SOS measure (2.2). For simplicity, suppose 
that the boundary height ^ at the right vertical side of A* is H' + 1 instead of 
= logL. There is an open {H' + l)-contour with endpoints at the endpoints 
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of the side with b.c. H' + 1. The probabihty that this contour equals 7 should 
be approximately given by the product of three factors: 

(i) the factor exp(— /3(|7| — L)) (the minimal length of the open contour is L 
and one pays for the excess length); 

(ii) a factor exp{+a{l3)A{^) / L) (with ^(7) the area to the right of the contour); 
this is due to the entropic repulsion and a{(3) should be approximately a(/3) = 
exp(— 4 X 2 X /3), where the factor 2 is due to the fact that H' + 1 — H = 2; 

(iii) exp(— 6(/3)^(7)/L) where b{(3) is approximately given (for /3 large) by 6(/3) = 
C2(/5) = exp(— 2/3) which appears in (6.1). 

Therefore, if /? is large the third term beats the second one and one pays both 
excess length and excess area, and it should be very unlikely that the contour 
reaches distance L j (log L) ^ \fL from the right rectangle side to which it is at- 
tached. We will find this probability to be roughly as small as exp(— cL/ (log -L)^), 
as would be the case for a random walk. Once we know the contour 7 does not 
go much farther than \/L away from the side of the rectangle, a suitable cou- 
pling argument will prove the theorem; see Section 7.1. Remark that without the 
Hamiltonian modification (6.1) (i.e. with fy = 0) the area gain kills the length 
penalization, and the contour would indeed invade the rectangle A*. 



7. Proof of Theorem 6.12 

The proof of Theorem 6.12 is based on the following lemma. Fix i = 1, . . . , M and set R := A*, 
R! := A*~^, so that the rectangle R\R' has horizontal length 2i, where i := L/(41ogL). Let 
also R" denote the rectangle of points in R at distance at least i from the right boundary. Note 
that Rd R" D R', and d{R \ R" , R') = i, see Figure 3. 

Let 7^(77), j = H' + 1, . . . , n+ denote the unique open j'-contour in the rectangle R enforced 
by the boundary conditions, attached to the right boundary. 



A, 





Ci 



Figure 3. Schematic drawing of the rectangles R, R" = R \ Aq, R' = R" \ Ai. 
Here the contour 7 = 'Jh'+i illustrates the event B in Lemma 7.1, while the 
chains Ci in Ai illustrate the vertical crossings used in the proof of Theorem 6.12. 
The shaded region corresponds to Int(7), while A(7) = i2\Int(7). The boundary 
between A2 and ^1 is dR' . 
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Lemma 7.1. Let B he the event that 7jif/_|_i(r/) does not intersect the rectangle R" . For every 
C > 0, n^/ [B") = 0{L-^) where n^'/ is as in Notation 2.1 and ^ is as in (6.32). 

We first show how to obtain Theorem 6.12 from the estimate in Lemma 7.1. The proof of 
Lemma 7.1 is given in Section 7.2. 

7.1. From Lemma 7.1 to Theorem 6.12. 

Reader's Guide 7.2. 

Let us first give a rough sketch of the couphng argument to be used. By condi- 
tioning on the value 7 of the contour ^h'+i one can roughly replace the measure 
TT^'/ appearing in Theorem 6.12 by the measure li^^^y where A(7) := A \ Int(7) 
is the region to the left of 7. Strictly speaking this is not true but we shall reduce 
to a similar situation by way of monotonicity arguments. Also, thanks to the 
argument of Lemma 3.2, one can neglect the influence of the ceiling constraint. 
Thus, one essentially wants to couple n^^^^ and 11^^ '■' on the region A*~ = R . 
Thanks to Lemma 7.1, one can assume tliat the rectangle Ai is contained in 
A(7). From the Markov property, it is sufficient to couple n^^^^ and 11^ on 
the interface separating the rectangles Ai and A2, see Figure 3. Thus, the de- 
sired estimate would follow if one could exhibit a coupling such that with large 
probability there exist chains €1,62 of sites in the rectangles Ai,A2 respectively, 
where both configurations are at constant height H'; see Figure 3. If there were 
no external fields and no wall constraint this would be a simple consequence of 
Lemma A. 2 (recall that, for /3 large, the interface is rigid and there is a density 
close to 1 of sites where the height equals the boundary height). However, due 
to the presence of the external fields and the floor at zero, establishing this fact 
requires extra work. The idea here is to reduce the effective size of the system by 
imposing boundary conditions H' on vertical crossings in the rectangle ^3, and 
in the rectangle Aq . More precisely, let p, pi denote two vertical crossings in ^3 , 
and let p2 denote a vertical crossing in Aq. Using monotonicity and the estimate 
of Lemma C.l of Appendix C, we replace H^^^^ by ^^'/^^) and 11^ by n^(p{^p2) ' 
where A(/3, 7) is the region between the chains p and 7, while A{pi,p2) denotes 
the region between the chain pi and the chain p2. Once this reduction has been 
achieved, the system is contained in the union of the four rectangles U^^qAj, a 
L X 4i rectangle, and one can easily show that since i is much smaller than L, 
and since H' = H + 1, the external field and the wall constraint can be neglected; 
see the proof of Lemma 7.3 below. At this point, one can use Lemma A. 2 to 
obtain the existence of chains €1,62 with the properties mentioned above. 

We turn to the details of the proof. It is sufficient to couple tt^'/ and vr^ on dR , the set of 
points in R' with a nearest neighbor in R\R' , i.e. 

IFa' ~ ^A Wr' = IFa' - ^A WdR'- 

Note that, because vr^'/ has maximal b.c. on a side of R and b.c. coinciding with that of 

TTj^ '■' on the other sides, vr^'/ stochastically dominates vr^^ on dR and therefore, by a union 
bound, one has 



Fa 
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where we define the events Ux,v '■= {ilx > v}. Next, we remove the ceiHng constraint from the 
measures vr^'/ ,7r^ . Since Ux,v are monotone events, we can estimate vr^'/ {Ux,v) ^ n^/ {Ux,v)- 

Moreover, as in Section 6.3, one has vr^ {Ux,v) = {Ux,v) + 0{L~'~^) where C is as large as 
we wish provided /? is sufficiently large. 

Let^ jj,'^ denote the marginal on R" of 11^/ conditioned to have ^h'+i{v) — 7- Let Int(7) 
denote all sites enclosed by the contour 7 and the right boundary of R, cf. Figure 3. Let A~ 
denote the set of sites x £ R \ Int(7) that have either a nearest neighbor in Int(7), or a site 
at distance 1/2 in Int(7) in either the south-west or north-east direction. Since conditioning 
on jH'+iiil) — 7 forces all sites in A~ to be at height r]x ^ H' (recall Definition 3.3 of an 
/i-contour), by monotonicity one has fP{Ux,v) ^ fJ-'^iUx^v) if Z^''' denotes the marginal on R" of 
n^'/ conditioned to have height exactly H' on all sites x G A~. Writing B for the event that 
'jH'+iiil) does not intersect the rectangle R" , one has, uniformly in x,v: 

Tr^J{Ux,v) - vrf '^(t/,,,) ^ L^^ + U^J {B') + max^T'([/,,,) - nf '^(^7,.,,). (7.2) 

Lemma 7.1 says that U^J (fi^) = 0{L~^), so that we are left with the upper bound on ^j^'^ {Ux,v) ~ 

U2 '■' {Ux,v) for 7 G ^. We now implement the system reduction mentioned in the sketch of the 
proof above. 

Let Ai, i = 0, 1, 2, 3, denote the L x £ rectangles in R depicted in Figure 3. Write Aj, for the 
external top boundary of Ai, i.e. the set of sites x ^ R such that x has a nearest neighbor on 
the top side of the rectangle A^. Similarly, write A'^ for the external bottom boundary of A^. 
Call £{Ai) the set of Z^-bonds e such that e has at least one endpoint in Ai and at most one 
endpoint in A* U A^. A vertical crossing in Ai is a connected set C C £{Ai) that connects Aj 
and A^; see Figure 3. 

Let J-"- (resp. denote the event that there exists a vertical crossing C in A3 such that 
rjx ^ H' for all X G C (resp. a crossing C in A^ and a crossing C' in ^0 such that rjx ^ H' , 
X G C U C). On the event one may consider the leftmost vertical crossing in A^ with the 
required property, where leftmost is defined according to lexicographic order. From the Markov 
property of , and using monotonicity, 

f^"'{Ux,v) ^ ^C<{T1_) + maxf,^'P{Ux,v), 

p 

where fi'^''^ stands for the measure fi'^ conditioned to have height H' on p, and p ranges over all 
possible vertical crossings in A3. Similarly, on the event J^^ denote pi (resp. P2) the rightmost 
(resp. leftmost) crossing in Aq (resp. ^3) and write 

K'^iUx,,) K'^iK)) Q'"'''{Ux,v) > min QP^'P^{Ux,.) - K'^iK) (7-3) 

Pl,P2 P1,P2 

where Qp^'P^ stands for the measure conditioned to have height H on pi,p2, and pi,p2 

range over all possible vertical crossings in ^0,^43. Altogether, 

p^{Ux,v) - K'^{Ux,v) ^ f^^i^l) + nf '-^(^f ) + max \p^^P{Ux,,) - Q'''^''{Ux,v)\ ■ (7.4) 

It follows from Lemma C.l of Appendix C that p^{J-l) and '■^(-F+) are 0{e~^^ Notice 
that fi'^'P (resp. Qp^'P^") are SOS measures with exactly H' b.c. around the domain whose 
boundary is determined by p (resp. pi) on the left and by A~ (resp. P2) on the right. Such 



In the sequel of the proof we introduce local notation for various conditional marginals of the measures 
n^'/, in order to keep formulas readable. 
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domain has (by construction) horizontal size of order I and vertical size L. To simplify the 
notation, we shall write Q for [i^'^ ^ qp\,pi ^ 

We now turn our attention to vertical crossings in the rectangles Ai,A2. Consider the in- 
dependent coupling ¥ fi"' , Q on Ai L) A2. Writing {r],r]') for the corresponding random vari- 
ables, let Ai, i = 1,2 denote the event that there exists a vertical crossing C in Ai such that 
Ve?? = Vgjy' = for all bonds e with both endpoints in C. Note that if C is a vertical crossing in 
Ai as above, then r]c = 7]'^ = H' , because of the boundary conditions equal to H' on the top and 
bottom boundary of Ai, i = 1,2. On the event ^1 n A2, one may consider the leftmost vertical 
crossing C2 in A2 and the rightmost vertical crossing Ci in Ai. From the Markov property of 
the Gibbs measures fi'^, Q and the fact that Q(-|r?Ci = i]C2 = H') ^''^(■kci = 1IC2 = H') have 
the same marginal on dR' (observe that dR' is just at the boundary between Ai and A2), one 
obtains that 

WiU,,,) - Q{U,,,)\ ^ FiAl) + P(^i). (7.5) 

We shall focus on the event A^, since the event A2 can be treated in the same way. To 
estimate P(^^), we use the fact (see e.g. [34, Lemma 11.21]) that non-existence of a vertical 
crossing in Ai implies the existence of a horizontal dual crossing in Ai. More precisely, let A[ 
denote the right side of Ai, i.e. the set of dual bonds e' such that e' crosses an edge of the 
form e = {x,y) with x ^ Ai and y £ R \ R". Similarly, let A^ denote the left side of Ai. We 
say that a dual bond e' is in Ai if e' crosses a bond e G £{Ai). Then, the event Af implies 
that there exists a connected set T> of dual bonds e' in Ai which connects the lines A^ and 
A^, and such that for every e' G P either Ve'?/ 7^ or Ve'V 7^ 0. Here we use the notation 
Ve'?? := Vgf/ if e' is the dual bond that crosses e. Moreover, for a given D as above, there must 
be a set F C P such that \V\ ^ |P|/2 and such that either Ey := {Ve'T] ^ for all e' G V} or 
Fy ■= {Ve'V 7^ for all e' G V}. Thus, using a union bound, one obtains 

nA)^Yl E if^^Ev) + Q{Fv)) (7.6) 

V VgV: 
\V\ > \V\/2 

where the first sum is over all connected sets of dual bonds P connecting A^ and A\ as above. 
We will need the following lemma. 

Lemma 7.3. There exist constants C, c,I5q > independent of /3 such that, for every set V of 
dual bonds in A with \V\ ^ 1/2, one has for all /3 ^ /3o 

max{/i^(^v.), Q{Fv)] ^ Ce-'^'^I^L (7.7) 

Let us conclude the proof of Theorem 6.12 assuming for a moment the validity of Lemma 7.3. 
From (7.6), summing over the possible (connected) sets V, and using \'D\ ^ I ^ , for all 
e>0, if/3^/3o: 



'C/3|T/| 



k^eV:\V\=k VCV: 
\V\^ k/2 



2V'=^^-/2 ^ 2C J] 6'=e-^/^'=/2 

k'^eV:\V\=k k^e 

^ ^ o(exp(-Li-")). (7.8) 

Since the constants implied in (7.8) are uniform in x, v and the choice of 7 G 6, the claim of 
Theorem 6.12 follows from (7.2) and (7.1). It remains to prove Lemma 7.3. This is where the 
reduction from fi'^ to fi'^'^ and 11^ to Qf^'f^ becomes important. 
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Proof of Lemma 7.3. We shall prove the bound concerning ^'^ = fi'^'P only, since the same proof 
works for Q = Qf^ '''^ . Consider the region Aq C i? delimited on the left by p and on the right by 
7. Since p is a vertical crossing in A3, one has A\ C Aq. A crucial fact is that |Ao| ^ AlLI. Let 
as usual vr^J denote the SOS measure on Aq with boundary condition H' outside of Aq, with no 
floor, no ceiling and no external fields. From Lemma A. 2 one has 7r^J(£'v') ^ e"'^'^'/^ for any 
V . Thus, it suffices to show that 

p?{Ey)^Ce'''T:Z(Ey), (7.9) 

for some constant C independent of /5. Note that the external fields contribute with the term 
0^2; Z^xsAo /a; ^ ^'^ Hamiltonian, and therefore, at the price of a factor e^^ we can 
remove all external fields in our measure [x^ . Then 

^kS"nx ^ 0, Va; G Ao) 
Next, from the FKG inequality, one has 

^ 0, Vx G Ao) ^ n *Ao(^- ^ 0) ^ n (1 -Ce-^^"""^) , (7.11) 

xSAo xeAo 

where we use the equilibrium estimate vr^J(??x < 0) = T^\J(r]x > H') ^ Ce~^^^'; see Proposi- 
tion 3.9. Since e~^^^' = e~^^ /L, one has 

J](l_Ce-^'3^'))^C^V^^^ (7.12) 

xeAo 

for a suitable constant Ci > 0. The desired conclusion follows from (7.10). This ends the 

proof. □ 

7.2. Proof of Lemma 7.1. 

Reader's Guide 7.4. 

Roughly speaking, the proof of Lemma 7.1 works as follows. There are — 
H' open contours attached to the right side of R (call it r) and let 7^, j S 
{H' + l,...,?!"*"} denote the j-contour. First one proves that the n"''-contour 
cannot reach distance say L/(logL)^ from r. For this, one lifts from H' to 
n"*" — 1 the b.c. around the three sides of R different from r (this is allowed by 
monotonicity) . This way, there is now a single open contour and the estimate 
follows from Proposition B.l. Next, we want to prove that 7„+_i cannot reach 
distance L/(logL)^ from 7„+, i.e. distance 2L/(logL)^ from r. Morally the proof 
works as for the previous case, except that now the b.c. n"*" at r is replaced by the 
b.c. — 1 at 7„+. The argument is then repeated iteratively and the statement 
of the lemma follows when j = H' + 1. In practice, there are many additional 
difficulties, which is why the proof is so much involved. The main obstacles are 
the following: 

(1) Proposition B.l cannot be applied directly, because it holds when both the 
floor constraint rj ^ and the field / are absent. However, Proposition 7.7 
will show that (morally) the field compensates the effect of the floor (which 
would tend to push the contours away from r); 

(2) Once 7^ is fixed, it is not true that the next contour (i.e. 7j-i) sees boundary 
conditions j — 1 in a new domain determined by 7^. The point is that, from 
definition of contours, we only know that the heights just to its left are at 
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most J — 1, not exactly j — 1. We will use monotonicity to be able to change 
to j — 1-b.c. 

(3) Applying Proposition B.l as outlined above to estimate the probability of 
large deviations of 7j-i given 'jj requires that the right boundary of the 
system (that is the configuration of 7^), where b.c. are j — 1, is not too wild. 
In practice, one needs it to be a path connecting top and bottom of the 
rectangle R, with transversal fluctuations at most of order say for some 
small e. We will apply the results of Appendix C to infer that, indeed, to the 
left of 7j and not far away from it there is a chain of sites, with transversal 
fluctuations of the required order, where heights are exactly j — 1. 

We use a sort of induction on the index of the open contours 7j, j € {H' + 1, . . . , n^}, where 
n"*" = logL. Let ^0 denote the rightmost L x i rectangle inside R as in Figure 3, and write 
^0 = ^^=H'+i^3 where Bj are non-overlapping L x io rectangles, ordered from left to right, 
such that £() = i/{n'^ — H') ~ L/(4(log L)^). Every rectangle Bj is further divided into two 
non-overlapping rectangles Bj, B'j, ordered from left to right, such that i?j is a L x £1 rectangle 
with ii = L\ for some (arbitrarily) small 5 > 0, and i?| = Bj\Bj \s a, L x £2 rectangle, with 
■^2 = -^0 — ^1 ~ ^0; see Figure 4. 



Bl 



r 



h l2 

Figure 4. The rectangles Sj, B'j, with the associated vertical crossings pj, for 
j = — 1 and j = n"*". 

Define vertical crossings in a rectangle as in Section 7.1. For j G {H' -|- 1, . . . ,n"'"}, consider 
the event Bj that there exists a vertical crossing Cj in Bj , such that ijx ^ j — ^ for all x ^ Cj. In 
particular, on Bh'+i, there exists a vertical crossing Ch'+i in ^0 with r]x ^ H' for all x S Ch'+i- 
Thus B D Bh'+i, and it will be sufficient to estimate from above the probability n^/(;B|^,_(_j^). 
Clearly, 

n+ 

ni'/(^^,+i)^ n«'/(^jne,+i), 

j=H'+l 
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where ^^++1 denotes the whole probability space. On the event Sj+i, let Cj+i denote the 
rightmost vertical crossing C in Bj^-^ such that r]x ^ j, x £ C. By conditioning on the event 
{Cj+i = Pj+i}, and using that the events Bj are increasing, one has 

n+ 

U^JiBj,,^,)^ maxM,^.,,(B,^), (7.13) 

where Pj+i ranges over all possible vertical crossings in -Bj+i (for j = , it is understood that 
coincides with the right boundary of R), and p-p^^^ stands for the SOS Gibbs measure on 
the region Aq) C R defined as the set of sites x G -R to the left of the crossing Pj+i, with 

• boundary condition r]^ = j for x G Pj+i and rjx = j — 1 on all other boundary sites. 
Note that a portion of the boundary height has been lifted from H' to j — 1 ^ H' . The 
advantage is that, this way, there is a unique open contour under the measure fJ-pj^^, 
rather than j — H' of them ; 

• floor constraint r]x ^ 0; 

• external field 

J E f^J-^-H- (7.14) 

seAjj.) 

where we recall that f^j = exp(— /3j)l^^ ^-ff+j) cf- (6.1). Note that the fields in (6.1) 
with index different from j — 1 — H have been removed. This is allowed since the 
function r] i— )• fx,a{il) is decreasing. The effect of the field fx,a is to depress the area of 
the {a + H + l)-open contour and, since there is just one open contour, we need only the 
term with a = j — 1 — H . 

Lemma 7.1 is then a consequence of (7.13) and the following claim. 

Claim 7.5. For j ^ H' + 1, uniformly in the vertical crossing Pj+i in -^j+i o-nd for every C > 0, 

Pp^^,{B^) = 0{L~^). (7.15) 

Let 7j denote the unique open j-contour for a configuration r] in the ensemble Ppj^i- By 
construction, 7j is to the left of Bj^^, and it may intersect the rectangle Bj or even Bj. Let 
Ej denote the event that intersects Bj. Conditionally on the event Ej, the contour stays 
to the right of Bj, and the estimate pp.^-^{B'j\Ej) = 0(exp(— L"^)) follows from (C.2), which 
is applicable since the shorter side of R is at most of length L. Thus the claim (and hence 
Lemma 7.1) follows once we prove the following lemma. 

Lemma 7.6. For j ^ H' + 1, uniformly in the vertical crossing Pj+i in Bj_^_^, and for all C > 0: 

Pp^^,{E,) = 0{L~^). (7.16) 

We will actually give an upper bound of order exp(— L^~^) for every e > 0. 

Proof of Lemma 7.6. For this proof, the crossing /Cj+i in ^j+i is fixed, and we simply write p 
instead of Ppj^i- Fix a contour F and consider the event -jj = F. Set A_|_ = Int(F) n A(j), and 
A_ = Aq) \ A_|., so that F is the set of dual bonds separating A_ and A_|_ within A^^^ (with 
Int(F) defined a few lines after (7.1)). For any F one may write 

p{jj = F) « e-^l^lz,- A_Zj,A+. (7.17) 

Here, Zj^\_ (resp. -^j,A+) is the partition function of the SOS model on A_ (resp. A^), with 
floor at height 0, field as in (7.14), b.c. j — 1 on (9A„ (resp. b.c. j on (9A+) and with the extra 
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constraint that t]x ^ j — ^ for all x G Ap (resp. rj^ ^ j for all x £ Ap), where Ap (resp. Ap) 
is the set of x G A_ either at distance 1 from A_|_ (resp. A_) or at distance \/2 from a vertex 
y S A+ (resp. y G A_) in the south west or north east direction. These constraints are imposed 
by the definition of j-contour; see Definition 3.3. 

Next, let (resp. -2^+) denote the partition function of the SOS model on A_ (resp. A_|_) 
with b.c. 0, no floor and no external fields, with the constraint that rjx ^0 for all x £ Ap (resp. 
r/a; ^ 0, X S Ap). Let w+ be the corresponding Gibbs measures. With these definitions, one 
rewrites (7.17) as 

Ki,=r) = ^e-mzo^_zl x (7.18) 

X a;_ (^e^ ^^eA_ ^v.^o.^^ - (j - 1)) uj+ (e^ ^"^*+ ^"^ ^ ; 77 > - 

where Z is the normalization and /C = Cj_i_H = exp(— /3(j — 1 — -ff)), see (6.1). 

We first observe that the very same arguments of Proposition 3.9 proves that for all x E A-t, 

^±(7,, ^j)xe-4«, (7.19) 

that is C^^e^^^-' ^ ^jJ±{'>]x ^ j) ^ C e~^^^ for some absolute constant C > 0. This is pos- 
sible thanks to the fact that even in the presence of the constraints on Ap the arguments of 
Lemma 3.7 can be used without modifications. Next, let vr stand for the infinite volume limit of 
the SOS measure with zero boundary condition. Proposition 3.9 implies that 7r(?/o ^ j) ^< e~'^^^ . 
Moreover, we observe that there exist constants c, to > such that for any x G A± at distance 
at least t > to from the boundary dA±, for any k: 

\uj±{7]x > k) - 7r(r/o > k)\ e"^*. (7.20) 

Let us prove (7.20) in the case x G A_. The case x G A+ is obtained with the same argument. 
Thanks to the exponential decay of correlations for the zero-b.c. SOS model at large /? (see 
[11]), (7.20) is equivalent to the statement obtained by replacing vr by vrj(_. Observe that by 
monotonicity uj-{r]x ^ k) ^ tt^ [rjx ^ k). Next, by the same argument of Lemma 3.7, the oj-- 

probability of a contour 7, is bounded above by e~^^'^\ and thus with probability at most e~'^^ 
there is no chain C of heights all greater or equal to zero in a shell of width t/2 around x, and 
at distance larger than t/2 from x. On the other hand, if E is the event that such a chain exists 
then by monotonicity and decay of correlations one has uj~{r]x ^ k;E) ^ vr^ {r]x ^ k) -\- e~^^. 
This proves (7.20). 

We turn to a rough estimate that allows one to rule out very long contours. Namely, if G 
denotes the event that ^ L-"^"*"^", then for all /3 large enough 

li{G) = l-0{e~^'*'°). (7.21) 

In what follows we may fix eq > as small as we wish. To prove (7.21), observe that from a 
trivial bound on the external fields and the FKG property for uj± one has 



.5E.eA. 1.. ^0 . ^ ^ _ _ ^ -Q ^ - {j- 1)). (7.22) 



a;eA_ 



From (7.19) 

u.{i^x > -{3- 1)) ^ (1 - ce-^^"^^-^^) ^ exp(-c7L), 
since j ^ H. Then (7.22) is bounded below by e~^^ for some C > 0. The same estimate holds 
for the last term in (7.18), and therefore one has 

/i(7, = r) ^ z^(r)e^^, (7.23) 
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for some constant C > 0, where v is the probabihty measure on contours V given by 

i/(r) oce-^l^lz^_Z^^. (7.24) 

Notice that v is the distribution of the unique open contour of the SOS measure on A_ U A+ 
with no floor constraint, with Dobrushin boundary conditions, namely with b.c. r]y = or ijy = 1 
depending on whether y has a nearest neighbor in A_ or in A+ respectively. It follows from 
(B.l) that i'(r) has the standard form 

where the decoration term ^' satisfies |^'(r)| ^ ce~'^|r|. The usual Peierls' argument shows that 

z.(|7| ^Li+^«) = 0(e-^'+^°), (7.25) 

and (7.21) follows. 

Thanks to (7.21), we can now restrict the summation in the normalization Z in (7.18) to 
contours F G G. Define 

:=-|A_|7r(??o>0), := |A_| 7r(r?o < -(j - 1)) , (7.26) 

$+ := -|A+| 7r(r/o ^ - 1) , := |A+| 7r(?7o < -j). (7.27) 

Proposition 7.7. There exists a < 1 such that for allT £ G one has the expansions (with error 
terms uniform inV £ G): 

uj- [e^ ^-e^- ^ ° ; r/ ^ - (j - 1)) = exp (/C| A_ |/L - - ^_ + 0(L°)) (7.28) 
UJ+ {ei ^^^^+ ^""^-^77^ - = exp ($+ - + 0(L")) . (7.29) 

Let us conclude the proof of Lemma 7.6 assuming for the moment the validity of Proposition 
7.7. First, observe that the functions in (7.26) and (7.27) satisfy, for some a < 1, uniformly in 
F G G: 

= ^(|A+| - |A_|) 7r(r?o > 0) , (7.30) 
^_ + ^+ = |A_| 7r(7?o ^ j) + |A+| 7r(r/o ^ i + 1). (7.31) 
From (7.28)-(7.31), setting |A(j)| = |A_| + |A+|, 4(/3) = 7r(r/o ^ k), and 4(/3) = ^'5fc(/3): 

u- {e^ ^^e^- ^ ° ; r? ^ - (j - 1)) u;+ (^e^ ^^S'^- ^"^ ^ r/ ^ - 

= exp (^(/C(l - 5i(/3)) - J,(/3)) + ^(/C(-l + 25i(/3)) + 5,(/3) - 5,+i(/3)) + 0(L-)) . 
Observe that 

Therefore, for large f3 one sees that 

-/C ^ /C(-l + 25i(/3)) + ~5j{p) - ~5j+i{P) ^ - /C/2. 
Since the term proportional to |Aq)| is independent of F, it plays no role in (7.18). Therefore, 

(l^\n^^ ^n^ray^ EreiJ.nG ^^(r) exp ( - ^|A+|) 

fi{Ej I G) ^ exp(0(L")) X — — , (7.32) 

EreG^(r)exp(-^|A+|) 
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where we recall that Ej is the event in (7.17), G the event in (7.21) and v the measure in (7.24). 
At this point an upper bound on fJ.{Ej \ G) follows from (7.32) by neglecting the negative ex- 
ponent in the numerator and using Jensen's inequality for the denominator. Using also (7.25), 
this gives 

^i{E,\G)^v{E,)e^v{^<\K\)) + 0{L^)). (7.33) 

It follows then from Proposition B.l that for every (3 sufficiently large, for all e > 0, if L is large 
enough: 

v{Ej) ^ exp(-Li-^). (7.34) 
Essentially, under u the contour T behaves like a random walk and the event Ej imposes a large 
deviation of order L/(logL)^ which is much larger than the typical diffusive fluctuation ^/L. 
Moreover, again from Proposition B.l one has 

K|A+|) = 0(i^+^). (7.35) 
Then (7. 33), (7. 34) and (7.35) end the proof of Lemma 7.6. □ 

Proof of Proposition 7. 7. Let us start with the lower bounds. Using first Jensen's inequality and 
then the FKG property for a;_ one has 

a;__ (e^ ^-'f^- ^ " ; ^ - (j - 1 



^ exp [- ''-(Vx ^ 1 7? ^ - {j- 1))]uj^{t] ^ - {j- 1)) 
r 

^ exp ^— |A_| - ^_ - ^_ , (7.36) 



where 



:= 7 ^ w-(??x > I ?7 ^ - (j - 1)) , 



L 



Similarly, 
where 



:= - ^ log(l - a;_(7?, < -(j - 1))). 
Je I''- ^ -1 ; ^ - j) ^ exp [^>+ - §+] , (7.37) 



a;eA+ 



Prom (7.20) one sees that both |^'_ — ^'_| and 1^*+ — \&+| are 0(L"), for some a < 1, uniformly 
in P G G. Therefore the lower bound in (7.28) follows once we establish that on G 

=0(L"), (7.38) 

for some a > 0. To prove (7.38) we use the following comparison estimate. Let us consider the 
case |$+ — By PKG one has uj+{i]x ^ — 1 1 r/ ^ ~ j) ^ i^^+iVx ^ ~ !)• C)n the other hand, 
whenever x G A„ is at distance at least , for some 5 > 0, from dA^, then we claim that 

w+(r/, ^ - 1) ^ tj+(r/, ^ - 1 h ^ - j) + 0(L"~') . (7.39) 
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These observations and (7.20) are sufficient to prove (7.38). In turn, (7.39) is a consequence of 
the technique developed below, cf. the comment after (7.55). 

To prove the upper bounds in (7.28) and (7.29), observe that from the FKG property of u}± 
one has 



El 

uj- I e 



rE.,A_i..,o.^^ -(j-1)) ^u;-(e"^^^^-'''^"")a;-( J] V^-O-Dj' 

a;eA_ 

Rewriting 

E.eA_ 1.. . ^ ^ I A„ l) c._ ( n (1 - ^-)) 



xeA- 

K - 



where (/^j; := 1 — e L^')a:>o^ and setting ■^j, = the bound (7.28) is then implied by 

( n (1 - ^ exp(-$_ + 0(L")), (7.40) 

c^- ( n (1 - ^-)) ^ exp(-^I/„ + 0(L°)). (7.41) 

xeA_ 

Similarly, the bound (7.29) is implied by 

w+ 1 ^.)) ^ exp($+ + 0(L°)), (7.42) 

u+y V^..)) ^ exp(-^+ + 0(L°)), (7.43) 

xeA+ 

with the notation (px := 1 — e"L"'"''^<°, and ■(/'a; = l^^<_j. Below, we establish (7.40)-(7.43) 
and (7.39). All these estimates can be achieved once one has an approximate factorization of 
the measure a;+ on a mesoscopic scale L", u G (0, ^). To illustrate this point, consider the 
expression (7.40), and suppose the product is confined to a square with side L", contained 
in A_. Then 

11(1-'^^))= E (-i)'^''^-(n'^^ 

k ) 



< exp(-^ a._(r?, >0) + O(l2(2--i))), (7.44) 

where we have separated the contributions of sets A with | A| ^ 1 and l^j ^2, and used the fact 
that ipx = ^l,,^>o + 0{L~^). In particular if one could factorize (7.40) into a product of (7.44) 
over all Qu C A_, then the desired bound would follow using also (7.20). 

To implement this idea we use the following geometric construction. Partition I? into squares 
P with side r = L" + 2L^ , where < 6 < u < ^ (we assume for simplicity that L", are both 
integers) . Consider squares Q of side L" centered inside the squares P in such a way that each 
square Q is surrounded within P by a shell of thickness , see Figure 5. Define the set S of 
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dual bonds associated to a non-zero height gradient, cf. Appendix A. The set S is decomposed 
into connected components (clusters) S. We call I{5) the collection of clusters S in S such that 
|5| ^ . Note that a cluster may have a nonempty interior. 

Consider the set of sites V C A_ defined as what remains after we remove from A_ all 
clusters S in I (6) together with their interior. A square Q C ^ is called good if the square 
P D Q has empty intersection with 9A„ UX{6); see Figure 5. We write G for the collection 
of good squares Q. The crucial observation is that if Q £ G, then there exists a circuit C of 



Figure 5. A drawing of the region A„. In the background the squares P (dashed 
lines). The clusters inside represent the set I{6) and the shaded squares are the 
set G of good squares Q. 

bonds of surrounding Q and contained in the square P D Q, such that r]c = 0. To see this, 
observe that there must be a circuit C of bonds surrounding Q such that gradients of ij along the 
circuit are 0, since otherwise there would be a path of dual bonds connecting Q with P^ with 
cross-gradients different from zero, and therefore a cluster S with size larger than intersecting 
P. Now, this implies that r/ is constant on C, and this constant must be zero, since otherwise Q 
would belong to the interior of a cluster S of size larger than because of the zero boundary 
condition on (9A_. 

Next, we estimate 1 — (fx ^ 1 for all x which do not belong to some Q £ G = G{I{S))- 
Therefore, summing over all possible realizations W of I{5): 

a;_( n ^ ^^-(X(5) = T^)L._(n ll{l-V.)\m = w) 

x&A- w QegxeQ 

^ ^ ^_ {I{6) = W)llsupn^c{ll(^-V^)) (7.45) 

w Qeg ^ xeQ 

where, for an arbitrary circuit C surrounding Q within the square P D Q and with a slight 
abuse of notation, we write vr^ for the SOS equilibrium measure on the interior of the circuit C 
with zero boundary conditions (without floor, ceiling and no fields). With the same argument 
of (7.44) one has, uniformly in C 

^c(n(l-^-)) ^ exp(-^ j;7r0(r?.>0) + O(L2(2"-i))). (7.46) 
xgQ xeQ 

Let Z{6) = be fixed. For any fixed square Q £ G, let Q' C Q he the square centered inside Q 
in such a way that Q' is surrounded by a shell of thickness within Q. Thus, if x £ Q' , then 
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X is at distance at least L from C, and therefore as in (7.20), for any p > 0, uniformly in C: 

4iVx > 0) = ffivx > 0) + 0{L~P) , xgQ'. (7.47) 

Prom (7.47), 

( n (1 - '^-)) ^ exp ( - ^ ^ 7r(r?,, > 0) + 0(L2(2«-i)) + ©(L-^'+^'-i)) . (7.48) 
There are at most 0(L^~^") squares Q. Therefore, from (7.45) one obtains 

^- ( n (1 - ^-)) ^ E ^- (^(^) = ^) ( - 1 E E *(^- > 0) + • (7-49) 

Next, we need to add back the contributions to the exponent in (7.49) from all removed vertices, 
where each vertex contributes at most 1/L. The contribution of a single removed shell P\Q' is 
0(L^+"~^), and they are at most 0(L^~^"), so that all removed shells give at most 0{L^'^^~'^) = 
0{L°^) for Q < 1 since u > 5. To estimate the contribution from all other removed sites, we 
observe that a site can be removed if it belongs to a square P that intersects either the boundary 
9A„ or the clusters of X((5), or if it belongs to the interior of a cluster of X((5). If A{X{5)) denotes 
the total number of sites in the interior of the clusters S G 2^(5), then these contribute at 
most K.L~^ X A{X{5)). Moreover, one has at most L^" x \X{8)\ sites that can be removed from 
intersections withX((5). These contribute at most /CL^"~-'^|X(5)|. Finally, one estimates roughly 
by L^"|9A_| = ©(L^^+^+^o) the number of sites removed from squares intersecting (9A„, since 
on the event G one has |i9A„| = 0{L^~^^°). Thus, the contribution from the boundary squares 
is 0(L2«+£o) = 0(L°), if 2u + £q< 1. Therefore, using /C ^ 1: 

7 E E > 0) ^ a>_ - L-^A{Z{5)) - l2"-1|X(5)| + 0(L"). (7.50) 

Thus, we have obtained 

^- ( n (1 - ^■-)) ^ exp ( - + 0(L")) a;_ [exp {L-^A{T{5)) + L^^-^\T{5)\)] . (7.51) 
xeA_ 

If {Si}'!^i denotes the collection of clusters of 1(6), with |X(J)| = \ Si\, then 

A{I{5)) ^\[j2\S,\y, (7.52) 

i i 

where the last bound follows from ^^xf ^ (X^j^O^ for all Xi ^ 0. Recalling that \Si\ ^ L'' 
for all i, using (A. 2), letting m represent the number of clusters 5i, . . . , Sm, and summing over 
their starting points xi, . . . , Xm, one has the estimate 

--(Ei^'^i^^)^ E E E ••• E 

i m ^ I Xl,...,Xm SlBXl Sm^Xm 
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where xi^i, ■ ■ ■ , Sm) = 1 if l^jl ^ for all i = l,...,m and IS*!! + •••|-S'm| ^ k, and 
x(<S'i, . . . , Sm) = otherwise. Therefore, 

i m^l \ X sbx, \S\ ^ 



m 



^ Ce-P^'^ E E C^^-^^'^ ^ e-^^'^ , (7.53) 

for any /3 large enough and for all L sufficiently large. From (7.52) and (7.53) one has 

i 

for all £ > and therefore 

(J- [exp (2L-M(X(5))] ^ ^ exp {21/ L - /5\/£/4) ^ + 1, (7.54) 

£=0 

since 2£/L ^ j3^/l./4:, for /5 large and £ ^ . Using (7.54), a Cauchy-Schwarz inequality and 
(7.53), it follows that 

a;_ [exp {L~^A{X{5)) + L^''-^\I{6)\)] ^CL^ e^(^") (7.55) 

for any a < 1. This ends the proof of (7.40). To prove (7.42) one repeats the same argument 
with the region A_ replaced by A+. 

We turn to the proof of the estimates (7.41) and (7.43). A minor modification of the same 
argument proves also the inequality (7.39). Here one has to replace the expansion (7.46) by the 
following bound: 

^c(n(l-^-)) ^ exp(- j;7r0(^,) + O(L-i+2-+'=(^)) + O(L6«-3))^ (7.55) 

where ^/^^ = 1^^ ^ j and c(/3) > can be made arbitrarily small by taking (3 large enough. Once 
this estimate (together with the corresponding statement for ip^ = Ir^^yj) is available, it is not 
hard to check that exactly the same arguments we used to prove (7.40) and (7.42) allow one to 
conclude. Here the term lCL~^7r{r]x > 0) appearing in (7.50) must be replaced by ■^{rix ^ j), 
which (thanks to j H + 1) is again less than L^^ for /3 large enough by (7.19). In particular, 
one can use the argument in (7.54)-(7.55) to conclude as above. 

It remains to prove (7.56). We cannot proceed as in (7.44) since ipx is not pointwise 0(1/L). 
From Bonferroni's inequality (inclusion-exclusion principle), one has 

4{ii(^-i^x)) = ^{-ir^4{u^^ 

xeQ Acq xeA 

^l-Yl ^"(V'x) + ^ E ^ci^-^y)- (7.57) 

xeQ x,y£Q: 
x^y 

Next, observe that 

E ^cii^^^v) = ( E ^c(^.))' + E + 0{L-'+'^) (7.58) 

x,yeQ: xeQ x,yeQ: 

x^y x^y 
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where Ti:^{'ipx]'^y) '■= '^ci'^^'^y) ~ '^ci'^^^^ci'^y)' ^^"^ ^c(V'x) = 0{1/L). We need the 

following bound. For some c(/3) — )■ as /3 — )• oo, one has 

x,yeQ: 

The bound (7.56) follows immediately from (7.57)-(7.59) and the fact that C^xeQ ^ci'^f^^))^ ~ 
0(L6«-3). 

To prove (7.59), first notice that by exponential decay of correlations [11]: 
for some constants ci, C2 > 0. Therefore (7.59) follows if we prove that for any constant C > 0, 

4{^,-i.y) = 0{L~^l^^<P^ 

0^|y| ClogL 

In particular, it suffices to show that uniformly in y ^ 0: 

7rO(V'o^,) = 0(^^'/'+^'(^)) (7.60) 

for some constant c'(/3) — )• as /3 — )• oo. The proof of (7.60) goes as follows. Let denote the 
event that there exists some fc-contour 7 that contains both 0, y. Then Ek-\-i C E)., 1, and 

j 

T^ci'^oipy) = T^ci'^oipy-^ Ef) + ^ TT^{il^oipy; E^ n E^^^) + 7r'](V'oV'?/; ^i+i)- 

k=l 

Now, if "ipQipy n El occurs, then there must be two separate families of nested contours reaching 
level J, one around and the other around y. By repeating the argument in the proof of 
Proposition 3.9, one has 

7:^c{i^o^y-E'i) = 0{e-^^^) = 0{L-^). 

If i/jQi/jyCi EkCi E^^-^ occurs then there must be nested contours around and around y separately 
from level A; + 1 to level j and there must be nested contours from level 1 to level k comprising 
both and y. In this case the argument in the proof of Proposition 3.9 yields 

7rO(V^o^,; Ek n i?,^+i) = 0(e-8/3(^'-^)e-^(i-^(^))^-^'') 

where iy denotes the length of the shortest contour comprising both 0, y and c(/3) decays as fast 
as 1//3. Since iy ^ 6, the above expression is 0{e~^^^^~'^^^'^^^) = 0{L~^^'^~^^'^^'^) for every k ^ j. 
Finally, the same argument shows that 

7:'c{^0i^y;E,+,) = 0(e-6/5(l-^(^))^') ^ q(^-3/2+c'(/3))_ 

These estimates imply (7.60). This ends the proof of Proposition 7.7. □ 

8. Mixing time in absence of entropic repulsion: Proof of Theorem 3 

Like for Theorem 1, it is sufficient to give the proof when n"*" = log L, the general case following 
easily (one needs to generalize the approach of Section 6.3 in the obvious way). For simplicity 
of notations, we call the SOS equilibrium measure with zero b.c. on dA^ and floor/ceiling at 
lb log L simply vr. Recall the definition (5.1) of the diagonal lines Ri and define, for 1 ^ j ^ N = 
(2L — 1)/L^/^"*"'' (assume for simplicity that N and L^^'^'^^ are integers), the subset Wj of as 
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Note that Wj D VFj+i is a roughly rectangular-shaped region of smaller side of order L^/^"*"^. 
Let also Sj denote the "brick" with horizontal projection Wj and floor/ceiling at it log L. From 
(2.8) and symmetry, we know that we have only to show that 

ll/i? - 7r|| ^ (8.1) 

with T = exp(c/3L-^/^+^^) and some large constant c. 

The proof is somewhat similar (but definitely simpler) to that of Lemma 6.6, so we will be 
very sketchy. The simplification is that, since the floor at — logL has essentially no effect at 
equilibrium, it is not necessary to introduce the field term (6.1) to compensate the entropic 
repulsion. 

We apply Theorem 2.2 with the following censoring protocol. We let AT = T/N and we let 
evolve first the brick Si for a time-lag AT, then S2 for another time-lag AT, and so on up to 
Sn. From Proposition 2.3 (which is immediately adapted to the case where A is not exactly a 
Lxm rectangle but rather is included in some, possibly tilted, Lxm rectangle) we have that the 
mixing time in each brick Sj, uniformly on the b.c. around it, is exp{0{(3n~^ L^/^"*"^)). Therefore, 
if c in the definition of T is sufficiently large, we can assume (modulo a negligible error term) 
that after the j*^ time-lag the j*^' brick is exactly at equilibrium, with b.c. on dWj n dAi, b.c. 
n+ = log L on dWj n Wj+i and, on dWj H Wj-i, a b.c. determined by the result of the evolution 
in the (j — 1)*'' time-lag. Theorem 2.2 then guarantees that the l.h.s. of (8.1) is smaller than 
ll/zy — vr||, with /z^- the law at time T of the censored dynamics. The inequality (8.1) then follows 
(via a repeated application of DLR) provided that one proves that, if ttj denotes the equilibrium 
on C/j- = TVi U • • • U Wj with b.c. on dUj HOAl and n+ b.c. on dUj D Wj+i, then 

llvTj - 7r\\u,., = 0{L-^) (8.2) 

i.e. the marginals of the two measures on Uj-i are very close. 

In analogy with the way Theorem 6.12 follows from Lemma 7.1 (cf. Section 7.1), to get (8.2) 
it is sufficient to prove that the open 1-contour does not intersect Wj-i, except with probability 
0{L~^). In turn (and in analogy to how Lemma 7.1 follows from Lemma 7.6), the desired upper 
bound on the deviation of the 1-contour follows if we prove the following. Consider a diagonal 
line Ri, with i ^ L^/^^^. Let A' = Ua ^ iRa and let p be a chain of sites in A', connecting two 
adjacent sides of A^;,, and at distance at most U from see Figure 6. 
















R_ 






Al 



Figure 6. A drawing of the triangular region A' and of the chain p. 
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The chain p disconnects Kl into two subsets and call A_ the one containing the North- West 
corner of A^. Let vr' be the SOS measure on A_, with b.c. on 9A_ n dA.^ and 1 b.c. on p. 
Then, the vr'-probability that the unique 1-contour reaches distance L^/^"*"^ from p is smaller 
than any inverse power of L. 

This is much easier to prove than the somewhat similar estimate of Lemma 7.6. The reason 
is that, since the fields (6.1) are absent and the floor has a negligible effect (recall that the floor 
was instead at height zero in Lemma 7.6), the desired estimate follows directly from a suitable 
modification of Proposition B.l, where the square Ql is replaced by a triangular domain. □ 

Appendix A. Peierls' estimates and low-temperature expansion 

Here we collect some rather standard facts concerning the low temperature expansion of the 
SOS model. With a small abuse of notation let Z\ be the partition function corresponding to 
the measure vr^. Following [11, Section 2] we will write Z\ as a sum over compatible cluster 
configurations. 

Definition A.l. A cluster X is a tuple (7, /ii, . . . , /i|^|), where ^ is a finite connected set of dual 
lattice bonds, and hi ^'L \ {0}. A cluster configuration is a collection of clusters {Xi, . . . , X^}. 

Let ^\ be the set of height functions r/ G with ?7x = for every x ^ A. Given 77 G 0^ 
one can define the associated cluster configuration {Xi, . . . , X^}, rn = m{r}), as follows. Fix an 
arbitrary orientation of the edges e = {x,y) of Z^. Let S = S{r]) be the collection of all dual 
edges e' such that the gradient of rj along the edge e = {x, y) crossing e' satisfies he ■= rjy — Tjx / 0. 
Let 71, ... , 7m denote the connected components of S. For each j = 1, . . . , m let Xj = (7^, {he}) 
denote the associated cluster, where {he} denotes the collection of gradients of r] along edges e 
that cross a dual edge e' G 7j . 

We define C{A) = U^^g^o {Xi, . . . , to be the collection of all possible clusters. Two 
clusters X,X' are called compatible, in symbols X ~ X' , iff 7 U 7' is not a connected set of 
dual edges, where 7, 7' denote the geometric part of X, X' respectively. Otherwise, X, X' are 
said to be incompatible, in symbols X r,^ X' . Also, let P(A) denote the collection of all pairwise 
compatible cluster configurations, i.e. of configurations {Xi, . . . ,Xm} with m ^ 0, G >C(A) 
for i = 1, . . . , m and Xi compatible with Xj for every i ^ j- Then, one has 

m 

E H-^^^i)' p(.X,)=exp[-f3j2\he\] (A.l) 

{Xi,...,Xm}eT'(A)j=l e 

where the sum over e extends over all {jjl edges e which cross a dual edge e' G jj. 

A.l. Peierls' estimate. As above, S denotes the random set of dual edges crossing a nonzero 
gradient. 

Lemma A. 2. There exists /3o > such that for all /3 ^ /3o, for all finite connected A C 1? , and 
all set V of dual edges, 

ttHSdV) ^e-(^-*)l^l. (A.2) 

Proof. We suppose that V is connected, since the general case follows by a standard general- 
ization. Let e' be a dual edge in V and let 5o denote the largest connected component of S 
containing e' . Then 

S:SDV 

where the sum is over all connected sets S of dual edges, such that S D V. Any r/ G such that 
So = S corresponds to a cluster configuration {Xs,Xi, . . . ,Xm} G 'C(A), where Xs is a cluster 
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of the form Xs = {S,hi, . . . ,h\s\)- For a fixed S one has T.hi^o,...,h^siy^o P(^s) ^ (4e if 
/3 ^ log 2. Therefore, using (A.l), neglecting the constraints on Xs, one has 

4(5o = S) ^(4e-^)l^l. (A.3) 

Summing over all S as above and estimating by the number of connected S 3 e' with jS"! = £ 
gives 

e^\v\ 

□ 

A. 2. Cluster expansion. We shall use a standard expansion for partition functions, adapted 
from [23,37]. For U a subset of the inner boundary of A, we write Z\^ij for the partition 
function with the sum over r] restricted to those r/ € such that t?x ^ for all x G U. 
One can write Z\^u similarly to (A.l): one defines C{A,U) as the set of all possible clusters 
(arising from height configurations respecting the positivity constraint in U) and P(A, U) as the 
collection of all pairwise compatible cluster configurations, with the same notion of compatibility 
as before. Then, one can check that (A.l) holds for Za,u, just with 'E'(A) replaced by V{A, U). 
We emphasize that it is here that one uses that f7 is a subset of the inner boundary of A: the 
identity would be false e.g. if U were the whole A. 

Lemma A.3. There exists /3o such that for all /3 ^ /3o, for all finite connected K (Zl? and any 
subset of its inner boundary U C A, 

logZA,u=J2^u{V), (A.4) 
VcA 

where the potentials ^ui^) satisfy 
(i) '^uiy) =OifVis not connected. 

(a) ipu{V) = (po{V) if dist{V, U) ^ 0, for some shift invariant potential V i— t- (po{V) that is 

My) = My+x) vxgz2. 

(Hi) There exists a constant (Sq > such that 

supsup|(^c/(y)| ^ eM-W-Po)d{V)) 
adv u 

where d{V) is the cardinality of the smallest connected set of bonds of 7? containing all the 
boundary bonds ofV (i.e., bonds connecting V to V^). 

Proof. We shall apply the main theorem from [.'i?]. Following [37], we define C as the set of all 
cluster configurations C that cannot be decomposed as C = Ci U C2 with two nonempty cluster 
configurations Ci,C2 such that {Xi,X2} is compatible for every Xi € Ci and X2 € C2. For a 
cluster X = (7, /ii, . . . , /i|-y|), define the function a{X) = \\X\, where A > is to be specified 

later and \X\ := Y^^li \hi\. Note that, for a fixed X = (7, hi, . . . , /i|^|), one has 

Y: e^^l^'I^X') ^ c(/3,A)e-(/'-2^)IVUc'(/3,A)|7|, (A.5) 

X': X'r^X 7'; 7U7'conncctcd 

where e.g. c(/3. A) = 2(1 - e-(/3-2A))-i ^/(^^ ^) ^ 3e-(^-2^)c(/3, A). So if /3 ^ 2A + 1, and A 
is larger than some absolute value Aq, (A. 5) implies 

Y e'^'''V(^0 a{X) . (A.6) 

X': X'rf^X 
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Equation (A. 6) corresponds to Eq. (1) in [-ji]. The main theorem there then ahows one to write 

log Zk,u= ^(^)' (A-^) 

C:CCC{A,U) 

for a function $ on cluster configurations satisfying $(C) = if C ^ C and 

|^(C)|e"(^) ^a(X), (A.8) 

for every cluster X, where a{C) := Y17=i ^i-^i) if = {Xi, . . . ,X„} and the notation C X 
indicates that Xi '/^ X for some Xi S C. The potentials ^ depend on U but for lightness of 
notations we keep this implicit. Taking X to be the elementary unit square cluster such that 
|X| = 4 in (A.8) one finds in particular that for every cluster configuration C one has 

4e^''(^\ (A.9) 

To write Z\^ij as in (A. 4) we follow [2.'), Section 3.9]. For any cluster configuration C € C, 
C = {Xi, . . . , Xn} with Xi = {-Ji, hi, . . . , write Cg for the geometric part of C, i.e. Cg = 

(71, . . . , 7„). For any G := (71, . . . , 7„) define 

ceC: Cg=G 

Using (A.9), if A ^ Ao, one has 

|V(G)| < 4e-^^-lT^I. (A.IO) 

Finally, set 

G=(7l,...,7n): 
Ujnt7i=\/ 

From (A. 7) one obtains the expansion (A. 4). The properties (i)-(ii)-(iii) follow as in [23] from 
an explicit representation of the function $(C), and from the exponential decay (A.IO). □ 

A. 3. Distribution of an open contour. Here we apply the expansion of Lemma A. 3 to derive 
an expression for the law of an open contour in the presence of a stepped boundary condition. 
Suppose a finite connected A C is given together with a boundary condition ^ with values 
in {0, 1} and such that it induces a unique open 1-contour 7. If 7 = F, for some connected 
set of dual edges F, then A is partitioned into two connected regions A+jA. separated by F. 
Moreover, 

4(7 = F) « e-^|riz^__^-Z^^^^H-, (A.12) 

where Ap are the sets defined after (7.17), and we use the notation Z\^i/ that was introduced 
in Lemma A. 3. By expanding the partition functions as in (A. 4), and retaining only terms 
depending on F, one finds that 

7ri(7 = F)ocexp(-/3|F|+^'A(F)), (A.13) 

where 

^a(F) = - Y Y ^At(^)+ Y V'A-(^)- 

VcA VcA+ VcA- 

ynr^0 ynr^0 vnr^fH 
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Here, the notation VCiT simply means that VCi (Ap U Ap) ^ 0. It is convenient to rewrite 
this expansion in the form 

^A(r)= Yl -^(^JT) (A.14) 

VcA 

where the "decorations" {4'{V;T)}vcA satisfy (cf. Lemma A. 3): 

(i) (/>(F; r) = if F is not connected. 

(ii) (p is shift invariant in the sense that 

(/.(F; r) = (/.(F + x; r + x) VxEZ^. 

(iii) There exists a constant /3o > such that 

sup|</>(y;r)| ^ exp(-(/3-/3o)d(y)) 
r 

where d{V) is defined as in Lemma A. 3. 

It is standard to check that these properties imply the existence of /3o such that, for any /3 ^ /3o 
and any i ^ I, 

V sup|,/.(y;r)| ^ exp(-(/3-/3oK). (A.15) 
— V 

d{v) > e 

Appendix B. Large deviations of the contour 

We begin by fixing some notation. S = {l,2...,L}xZ will denote the infinite vertical strip of 
width L. We denote by A, B the points of coordinates (1, 0) and (L, 0) respectively. The L x L 
square with corners A, B, C, D, where C = (L, L) and D = (1, L) will be denoted by Ql. Next 
we fix an open contour L^ inside Ql joining A with B with the property that stays above the 
line at zero height and does not reach height L^, for some 6 < 1/2 that in the applications will 
be taken small. The region inside S above is denoted by A and we set Q = n A. We let i^q 
be the law of the open 1-contour T joining A with B, for the SOS model without floor/ceiling 
in Q, with 1 b.c. along and b.c. otherwise. We know that vq can be written as 

i/q(F) « exp(-/3|F| + >I/q(F)) (B.l) 

where is the function appearing in (A.14). Fix a G (1/2,1) and i € [L", L/ log(L)^] and 
define Eg as the event that the path F reaches height £ (note that £ ^ L^). 

Proposition B.l. Uniformly in as above, there exists /3o independent of {£, L) such that, 
for all 13 > (3^ and all L large enough 

UQ{Et) ^ c' eM-c£^/L) 

for some constants c,c' . 

B.l. Proof of Proposition B.l. As a first preliminary step we remove the dependence on the 
upper boundary of Q. Let u\ be the probability distribution on contours in A joining A, B given 

by 

i/a(F) ocexp(-/3|F| + ^A(F)). 

Claim B.2. For any (3 large enough 

UQ{Ee) ^ 3uA{Ei) + e-''^ 

for a suitable constant c = c(/3). 
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Proof of the Claim. Let Qq and Qi be the set of contours which stay below height L — log(L)^ 
and height L respectively. Then 

"^^^^^ ^ ^A(x.oe^-(^)) + .A(x.oe^*(^)) ^ 

where 

A^(r) = ^Q(r) - ^A(r) 

and the inequality sign comes from restricting the average in the denominator from contours in 
Gi to contours Gq. Since minregg |r| — minre^o l-'^l ^ L, a, standard Peierls argument shows that 

for some /3o. Moreover, thanks to the exponential decay of the decorations (A. 14), 
|A^(r)| ^ i VrG^o and |A^'(r)| ^ e'^'^-'^o^L ^reGi- 

Therefore the first term in the r.h.s. of (B.2) is smaller than 'iv\{E() while the second one is 
bounded from above by e~'^^ for some constant c = c(/3) diverging as /3 — oo. □ 

Back to the proof of the proposition: since the event is increasing, we can change the b.c. 
from to 1 along the lateral sides of A, up to height (3/4)^ (note that in this situation the 
endpoints of V are shifted upward by (3/4)£). We still call u\ the measure of F in this situation. 
Again by FKG, we have 

where is the increasing event that F stays at distance at least from F*, for some small 
but positive constant e. 

Thanks to the decay properties of the potentials (j){V]T), for every F in G"*" we can replace 
'I'q(F) with ^'§(F), up to a negligible error term. Then, the ratio (B.3) equals 

+ (1 + 0(1))^ (B.4) 

with the measure of the contour for SOS in the strip S, with b.c. above A + (0, (3/4)^), B + 
(0, (3/4)^) and 1 b.c. below it. 

Note that, if the complementary event [G~^Y happens, it means that the contour F makes 
a downward deviation at least (1/2)^ from its natural height (3/4)^. Therefore, f§(G+) ^ 1 — 
f§ (£"£). As a consequence, it suffices to prove: 

Claim B.3. For any c > and all (3 large enough depending on c, one has f§(£'^) ^ e~^^^/^. 

Proof of the Claim. By translation invariance, we can assume that the path F starts at A = 
(1, 0), ends at i? = (L, 0) and and replace Eg with Eg/^. We would like to appeal to the results 
of Section 4.15 in [23]. For this purpose we need to tackle the fact that decorations touching the 
boundary of S may behave differently from decorations inside S. We therefore introduce a third 
(!) probability measure on all paths F between A and B (even those going outside S) denoted 
simply by P(-) and corresponding to the weig ht e-^|r|+*z2(r) and we write 

f§(F reaches height £/4) 

e(f reaches height ^/4; F G S; e*s(r)-*z2(r)' 



E(FgS; e*s(r)-*z2(r) 
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Using Section 4.15 of [23] we get that 

P(r reaches height ^/4) ^ e""^^/-^ 
for some constant c > 0. On the other hand, (A. 15) implies that 

|^'22(r) - ^§(r)| ^ e-^^-^'''^\h e r : dist(6,S'=) ^ log(L)2| 

which implies 

]E('e2|*z2(r)-*s(r)|\ < gc'iog(L)2 



for some constant c', thanks to the large deviation results of Section 4.15 in [2.'!]. Finally, thanks 
to (A. 15) and Proposition 4.18 in [2.')], if C is the cigar-shaped region with tips at A,B defined 
by 

^ cp(r G c) ^ g-c"iog(L)2/-^ 

where the constant C is a deterministic lower bound on e*^*-'"^"*^^^^^ for F £ C obtained again 
using (A. 15). □ 

Appendix C. 

Fix a £ (0, 1). Let R be the intersection between and a, L x L"" rectangle, not necessarily 
parallel to the coordinate axes. Let A C be such that A contains R and is contained in some 
2L X 2L square. A subset C = {xi, X2, . . . , x^} of will be called a spanning chain if 

(i) Xj+i) = 1 for ah i = 1, ... A; — 1; 

(ii) C connects the two shorter sides of R. 

For a fixed n ^ let J-+ ( ) be the event that there exists a spanning chain where the surface 
height is at least (at most) n. 

Lemma C.l. For (5 large enough 

nx(n) ^ n:('^(^;) ^ e-^^^ (C.l) 

Assume moreover that £{A)e~^^^"'~^^^ ^ 1, where £{A) is the shortest side of the smallest rectangle 
containing A. Then 

Here, as in (6.1), the field is / = J2yeA fy 



n;^'^ (J-f ) ^ UliFl) s; e-"^". (C.2) 



Proof. We first observe that J^+(J^_) is an increasing (decreasing) event and therefore the first 
inequalities in (C.l), (C.2) are trivial because the fields fy are decreasing functions. Again by 
monotonicity H^-^ {J-"^) is bounded from above by the probability w.r.t. the SOS model tt^"^ 
without fioor. Moreover F^{J-'t) occurs iff there exists a *-chain {yi, . . . , y„} connecting the two 
long opposite sides of R and such that rjy. ^ n — 1 {r]y^ ^ n + 1) for aU z. In turn that imphes 
the existence of a (n — l)-contour ((n + l)-contour) larger than L". 
As in the proof of Lemma 3.7 we get that 

n'^^i'y is a (n - l)-contour) ^ (,-f^\^\+i^.eA-, IIMU ^ ^-§h\ 
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where in the last inequahty we used ||/x||oo ^ together with |A^| ^ 2L|7|. Simple counting 
of 7 finishes the proof of (C.l). 

Similarly, it follows from Proposition 3.6 that 

Uli^ is a (n + l)-contour) ^ e-'3|^l+^^-'''"+''l^-l 

Isoperimetry gives |A^| ^ £(A)|7| which, combined with the assumption l(A)e~'^^^"'~^^^ ^ 1, 
implies 

n^(7 is a (n + l)-contour) ^ e'^ltl 
and the proof of (C.2) follows. □ 

Appendix D. Proof of inequalities (6.22) and (6.29) 

Proof of Lemma 6.11. Fix i > 2. By removing the field / of (6.1) we only increase the 
surface so to bound the probability of the decreasing event we may work in the model vr^^ , 
i.e. the standard SOS model on A^, with no field and floor/ceiling at height O/n"*" = logL. If 
fails then for some R ^ TZ we can find contours {(7^, /is)}se,y satisfying the hypothesis of 
Proposition 3.6 each with |A^^ n i?| ^ 1 and hg H' + 1 such that {^s<^y^'ys,h holds, and that 

^ |A^^ n i?| ^ £L/2. (D.l) 
For a given ensemble of contours as above, define a sequence of subsets Wi C A by 

v^^o = A, 

W^= IJ A^^ for i = 1, . . . , n+ - fl-' . 

si^y ■.hs=H'+i 

Let A denote the set of all possible such collections of contours {(7^, /is)}. For all i ^ 0, let 

Let A{a) = ^(01,02, • ■ • ,an+-_ff') denote all collections of contours matching a given sequence 
of Oj's. Then (D.l) is equivalent to Ylii > 1 ^ (.L/2. Since i? is a diagonal, |A^^ n i?| ^ \\ls\ 
and so 

n+-H' 

\ls\ ^4 a,. (D.2) 

For any C A let 

B{W) = Y 

(7l,7^,...,7,'n) 

where the sum is over all collections of edge-disjoint contours {7I}, with pairwise disjoint interiors 
{A^/} all contained in W and with |A^/ D R\ ^ 1 for all i. Any such contour must have an edge 
adjacent to some v £ W n R in the dual lattice Z^*. If e is an edge in the dual lattice Z^* then 
there are at most 3" contours 7 of length n containing e. Hence for large enough /3, 

/ 00 \ mr\R\ 

BiW) < + S^e-^'^/^)" ^ exp(|I^ D R\) , 

since each contour must contain at least one edge adjacent to some v G WDRin the dual lattice 
Z^*, there are at most 4|VF n i?| such edges and the contours are edge-disjoint. 
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Now for {{js , hs)} s£.y S A{a) by Proposition 3.6 we have that 

^ exp(-|/3^|7,| J (D.3) 

for any (3 ^ Cq since e~^^^'' ^ ^ ^-i and |A^J ^ (L/4)|7s| for any contour 7^ by the 

isoperimetric inequality in T?. Substituting this expression we have that 



^ exp f - 2/3 ^ J ( - (^/^) E 



where the last inequality is by (D.2). This in turn is at most 

, n+-H' ^ n+-H' , n+-H' n+-H' ^ 

exp(-2/3 Y «0 n S (W^-i) ^ exp ( - 2/3 ^ + ^ a,„ij 

^ i=l ^ i=\ ^ j=l i=l ^ 

^ exp . 

The final inequality follows for large /3 since uq = L. As there are at most ^ 
choices for a = (ai, 02, . . . , an+-H') have that 

<(Gi)^i-E E <(n^7... 

^ 1-L'°g^exp (^-^/3-^i> 
^ 1-exp T-^^L 



for large /3, as required. □ 
Equation (6.22) follows similarly with a simpler proof. 
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